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ABSTRACT 


The kinetics of the oxidation of ferrocyanide by 
both horseradish peroxidase compounds I and II (HRP-I, 
HRP-~T TT) have beens studied: asia function. of pH at, 25°.»..The 


, for the reac- 
2app 
tion of HRP-I with ferrocyanide varied from 4.1 x ye 


ee es oe ee the BH xange, a. oO il. 3’. 


apparent second order rate constant, k 


The value of ee for the reaction of HRP-II with ferro- 
cyanide, varied. from. 1. 3x ters coe ten ek OCs Mee ce oie 
Over ice pi range 2). Oo. tO. 10.3... Whe; pH.dependence,of the 
HRP-I-ferrocyanide reaction has been interpreted in terms 
ofa ee ae Ses Gi ch oso) at tne ace Va site of 
feo we Mow Cependence Of the HRP=IT—-ferrocyanide reac- 
EPon has been interpreted in terms of three ionizations with 
DiS Of 5.4, ).27an0 826 at the active site of HRP-II. 

The steady state kinetics of the HRP catalyzed 


oxidation of ferrocyanide by H were studied by the 


203 
measurement of initial reaction velocities. Rate constants 
obtained by the steady state method were in agreement with 
those obtained from studies of isolated reactions on the 
stopped flow apparatus. The kinetic results can be correlated 
with conclusions drawn from previous kinetic studies on 

native HRP and with current ideas on the structures of HRP- 

I and HRP-I1. 


The kinetics of the binding of imidazole to ferri- 


protoporphyrin IX (hemin) in aqueous ethanol (44.5 weight %) 
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were studied at 25° using the temperature jump relaxation 
technique. The reaction was studied both as a function of 
imidazole concentration and hydrogen ion concentration. 

The acid dissociation constant for imidazolium ion was 
determined by acid-base titration and that of the protonated 
hemin species by spectrophotometric titration. The kinetic 
results may be explained quantitatively by a reaction in which 
imidazolium and imidazole bind at different rates to the 
basic hemin species. A comparison of the rates of ethanol 
exchange and imidazole binding to hemin favoured imidazolium 
ion Epes on on the hydroxy hemin species by a Sxl ion 


pair mechanism. 
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CHAPTER 1 
TnHtroduction 


1.01 Horseradish Peroxidase 

Horseradish peroxidase (EC 1.11.1.7; donor: 

H,05 oxidoreductase) is an enzyme that catalyzes the oxida- 
tion of certain substrates by methyl hydrogen peroxide, 
ethyl hydrogen peroxide and hydrogen peroxide. 

As early as 1810 the catalytic properties of fresh 
horseradish root were recognized (Planche, 1810). However, 
a proportional relationship between peroxidase activity 
and absorption in the Soret region was not shown until 1931 
Weuimeoe ale, t75l). Kedsorably pure preparations of *horse= 
radish peroxidase, HRP, were first obtained by Keilin and 
Bann tro37).. Peroridases are found in plants, animals and 
microorganisms. Horseradish roots, Japanese radish, turnips, 
enadetfg cap all Contain relatively fiigh coneentrations of 
peroxidase and serve as sources. Peroxidases are also 
present in milk, the thyroid gland, human saliva and blood 


plasma (Saunders et al., 1964). 


Purification and Separation 

A number of isozymes of HRP from horseradish roots 
have been observed by several workers (Theorell, 1942; Paul, 
£958; Kasinky and Hackett; 1967) Shannon et ‘al., 1966). The 
amino acid composition of the most abundant isozymes of HRP 


as determined by Shannon et al. (1966) are listed in Table 
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(1.01). All isozymes found for HRP contained ferriproto- 
porphyrin IX as their prosthetic group as do most of the 
other plant peroxidases. The molecular weight of HRP is 
approximately 40,000 (Keilin and Hartree, 1951). HRP is 
unusual as an enzyme because it is a glycoprotein. Carbo- 
hydrate analysis has shown that neutral and amino sugars 
account for 18% of the enzyme (Shannon et al., 1966). The 
neutral sugars in the different isozymes include galactose, 
arabinose, xylose, fucose and mannose. The amino sugars 
include mannosamine and galactosamine. The function or the 
mode of attachment of these sugars to the protein is as yet 
unknown. Separation. of the carbohydrate and protein portions 
of the molecule was not achieved throughout normal separation 
procedures (Paul, 1963). Keilin and Hartree (1951) described 
properties of HRP in its alkaline and reduced forms as well 
as the compounds it forms with carbon monoxide, cyanide, 
azide, fluoride and hydrogen peroxide. They also concluded 
that a water molecule normally occupies the sixth co-ordination 
POsSition.of the Fe(II1) ion,.of HRP. 

The purification procedure of commercial samples 


of HRP is described in Appendix 3. 


Reaction of HRP with Oxidizing Agents 

HRP was first observed to form spectroscopically 
distinct compounds with hydrogen peroxide over 30 years 
ago (Keilin and Mann, 1937). When H,0. is added to a dilute 
(107?) solution of a highly purified HRP preparation the 


brown colour initially present gives way to a green colour. 
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Table 1.01 


: a 
Amino Acid Composition of the HRP Isozymes 


Amino acid Isozyme C Isozyme B 
(gm residue/100 gm protein) 


lysine oI 1.8 
histidine 1.0 aa 
arginine Ciara on 
aspartic acid 14.5 IMS als: 
threonine De 6.6 
serine 4.7 Se! 
glutamic acid GieZ 6.6 
proline 4.0 ae 
glycine PRE ES 
alanine 4.0 4.4 
half-cystine ORS. eu 
isoleucine 2 a0 SI 
leucine EQ ak £0 
tyrosine AY te 
phenylalanine Diers 8.8 





2 from Shannon et al. (1966) 
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Upon standing for several minutes or on adding a small 
amount of reducing agent Such as p-cresol. the,green colour 
gives way to a red one. With an enzyme preparation that is 
less pure, the reducing agent present in the preparation 
(endogenous donor) causes the green compound to change 
immediately to the red one. When an excess of the reducing 
agent is added, the original brown colour of the solution is 
returned. The green compound that precedes the red one is 
known as horseradish peroxidase compound I and is abbreviated 
to HRP-t, HRP-I was first observed by Theorell (1941) who 
used HRP of reasonably high purity obtained by repeated 
fractional crystallization. The red compound following 
HRP-I formation is known as compound II or HRP-II. At 
higher hydrogen peroxide concentrations a bright red com- 
pound, HRP-III, is formed. Changes from HRP-I, HRP-II, and 
HRP-III to HRP are reversible as HRP can be reformed in the 
presence of an oxidizable substrate. Under conditions 
where HRP-III is present for a number of minutes the enzyme 
activity is reduced (Chance, 1952). At excessively high 
concentrations of hydrogen peroxide an irreversible change 
takes place with the formation of the emerald green com- 
pound IV, HRP-IV, (Chance, 1949a). Chance (1949b) showed 
that HRP-I preceeds HRP-II and that HRP-II is formed only 
from HRP-I... George (1952, 1953a) titrated HRP-IL with 
ferrocyanide and showed that HRP-II was converted to HRP 

by a one-electron reduction. Similarily Chance (1952b) 


showed that the conversion from HRP-I to HRP-II also involved 
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a one-electron reduction. Thus HRP-I retains the two oxidiz- 
ing equivalents of hydrogen peroxide and HRP-II only one. 
Only HRP=I and HRP-ITI are Siar to be enzymatically 
important intermediates in the peroxidatic oxidation of 
substrates. The relationship between each of these forms of 
the enzyme, illustrated with ferrocyanide as the substrate, 
are shown in Figure (1.01). The features of the absorption 
spectra of HRP, HRP-I, HRP-II, and HRP-III are listed in 
Table .¢&02). 

HRP-I may be formed from hydrogen peroxide, methyl 
hydrogen peroxide, and ethyl hydrogen peroxide (Chance, 1949b). 


The reaction of HRP and peroxide is first order in each of 


Bnew teacconts. For reaction (1.01) (Chance, 1963): 
k 
HRP + ROOH pee epi (101) 


where R = HOOH, MeOOH, EtOOH 


the second order rate constants for each of hydrogen peroxide, 
methyl hydrogen peroxide and ethyl hydrogen peroxide are: 

6 6 on eee =—1 ge 
pee Oe ele ae LO me Ox LO Mi sec  .. In, addition to whe 
three peroxides known to give rise to HRP-I and HRP-II, a 


number of strong oxbaiaing agents (HOCL, HOBr, Nacio Cr 


oe 2 


3! KIO, and 03) cause similar reactions (George, 1953b). 


KBrO 
Reaction with Reducing Agents 

The HRP - hydrogen peroxide system is known to 
oxidize inorganic ions such as ferrous ion, ferrocyanide, 


chloroiridete, nitrite, and todide (Paul, 1963) as well 


as numerous organic substrates including aniline, mesidine, 
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Table 1.02 


Spectroscopic Features of HRP and its Hydrogen Peroxide 


Compounds? ’??¢ 
HRP 
r 403 497 530 580 641 
max 
€ 91 ae eet feds) aS Zac 
mM 
RRE<t 
r 410 
max 
EM 48 
HRP-Il 
r 418 Lat Sens: 
max 
: Tig) 
eM 25 S35 
HRP-IITL 
r 416 546 583 
max 
EM 106 10 Ceo 


abe: md oe: 
With ee the band position in mu and EM the millimolar 
absorptivities 

2 From Paul. (1963) and Brill and Williams. (1961). 

eit addition, the following isosbestic points exist between 
each of the compounds (Chance, 1963a): 
Compounds Wavelength 
Hee 1, -onRP= LT 395 mu 
HRP - HRP-II 410 my 


HRP - HRP-I 427 my 
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p-cresol;, and tyrosine (Saunders -et al., 1964). 

Chance (1963a) studied the kinetics of the reaction 
of HRP-I and HRP-II with a number of oxidizable substrates 
Including p—aminobenzolc acid, nitrite, and ascorbic acid. 

In general the rate constant for the reaction of HRP-I 

with substrate is generally 40 to 100 times faster than for 
HRP-II. Thus the reaction of HRP-II with substrate is 
usually the overall rate determining step in the enzyme 
Feaccion. 9in the Case Of p-aminobenzoic acid acting as 
SUbstrarce, Sy and giving product, P, the rate constants for 


the individual reactions determined by Chance (1963a) are: 


HRD HO SOTRPSt k =notso 0b oscar (1762) 
eee lapp 
ARPS tech S th RDA DLE: IPIOtk . =2x 100 Moee oi (e038 
app 
e eed =) 
PRP 4, S > GarRpre! p sa = 2.4 x 10°M “sec (1.04) 


The second order rate constants are subscripted as apparent 
rate constants as they may show considerable pH dependence. 

A recent stopped-flow study of the hydrogen 
peroxide oxidation of luminol by HRP has indicated the second 


order apparent rate constants, k and k EO be. 2a 


2app 3app’ 


TO ee ens, ayia ade O (Cormier and Prichard, 


1968) at pH 8.0. 


Reaction of HRP with Ligands 

HRP reacts with a variety of ligands including 
cyanide, fluoride, azide, hydrogen sulfide, nitric oxide, 
and hydroxide (Keilin and Hartree, 1951). The kinetics 


of ligand binding to HRP has been used as a probe to gain 
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information about proton ionizations at the active site of 
the enzyme (Dunford and Alberty, 1967; Ellis and Dunford, 


196.3)". 


Ionizable Groups at the Active Site of HRP 

A number of studies have measured proton ionizations 
associated with the active site of HRP. Values of the acid 
dissociation constants are listed as pK's in Table (1.03). 
The physical techniques used to measure these ionization 
constants include spectrophotometric, magnetometric, and 
potentiometric titration as well as kinetic methods. It 
is generally agreed that the pK values on HRP in the range 
10.8 to 11.3 are-associated with the formation of an hy- 
Pea complex on the single iron atom of HRP (Elizvs and 
Dunford, 1969). No general agreement exists as to’ the 


nature or even the presence of the other ionizations. 


1.02 Hemin 

The porphyrin ring system has unique 
chelating properties towards eats metal cations, 
The metalloporphyrins constitute the functional group 
of a wide range of biologically important compounds. 
Ferrrprotoporphyrin LX or hemin; an i1ron(lit) porphyrin, 
which is the prosthetic group of HRP, is shown in Figure 
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The pK Values Found on Various Species of HRP 


Species pk (2) u Method 
HRP® a5, 0 = x: Spec 

HRP Omg t 20 noc constant Spec 

HRP T27 20 ~0.06 Mag 

HRP 10.56 30 0.07-0.19 Pot 

HRPo 11.0 25.0 Antal Spec 
HRP Bete BA a1 25.0 0.11 Kin-CN 
HRpo't 7 iade et at ctean 25.0 0.11 Kin-F" 
HRP-CN© 6. Tail 0.87 25.0 Orit Kin-CN_ 
HRP-F° 5.2 25.0 (ieee: Kin-F 
HRP-112 3.4, 5.2, 48.6 25.0 Omid Kin-Fe (CN) ¢ 
HRP-19 5.3 25.0 Oe it Kin-Fe (CN) -~ 
Spec: spectrophotometric equilibrium titration 

Mag: Magnetometric equilibrium titration 

EOC: potentiometric equilibrium Poe eon 

Kin-CN : kinetically, cyanide binding. 

Kin-F i kinetically, fluoride binding 

Kin-Fe (CN) ¢: kinetically, Fe (CN) - substrate 

= (Theorell and Paul, 1944); b (Hneore ll 7 1942) 3 (Hac DULY, L057) 
d (ellis and Dunford, 1969); © (Ellis and Dunford, 1968); 


(Dunford and Alberty, 1967); J see Chavter) 35 
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Figure 1.02: Ferriprotoporphyrin IX or Hemin. 
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SELUCCUre 
The chelation of a metal ion by a porphyrin 
involves the displacement of two protons from the nitrogen 


atoms to form a structure that is almost square planar. 


ae Moe 
Figure 1.03: Placement of the Nitrogen Atoms About the 


Central Metal Ion. 


The fifth and sixth co-ordination positions above and 

below the plane of the nitrogen atoms are available mae 
further co-ordination to form octahedral or square pyrimidal 
structures. ‘The X-ray structural determination of chloro- 
hemin, the mono chloro derivative of ferriprotoporphyrin IX, 
has indicated that the compound is nearly square pyrimidal 
with a 0.48 A out-of-plane displacement of the iron atom 


°o 
(Koenig, 1965). The Fe-Cl distance was found to be 2.22 A 


12 


= : a a 
; *) yee 


a 
Pn eet: ec | 
& tT SoOU iso 
a ne 


‘Tad y S161 er *O. 


. 





La 


and the Fe-N distance 2.00 A. Le BZ=ray structure of “the 
bis (imidazole) -a,8,y,é-tetraphenyl-porphinato iron(IIT) 
chloride has been determined (Countryman Strele; "1709) “alia 
the environment of the iron atom is more nearly octahedral 
than the chlorohemin compound with only a slight (0.009 A) 
out-of-plane displacement. Magnetic susceptibility 
measurements Of chlorohemin (Pauling and Coryell, 1936) 

and aquohydroxyferriproto-porphyrin IX (Coryell et al. 
1937) have indicated that these compounds are mainly 
Peesenitein their high spin state. A number of other solid 
hemins with ligands such as bromide, formate, acetate, azide 
and propionate are also high spin (Rawlinson and Scutt, 
1952). The solid bis(cyanide) and the bis (imidazole) hemins 


ere) Ow spin complexes {Countryman et al., 1969). 


Hemin as the Prosthetic Group of Heme Proteins 
Thewterriceand erro protoporphyrin IX's are 
found in a number of heme proteins (Falk, 1964) and among 
these are the hemoglobins, ‘myoglobins, peroxidases, cata- 
lases, and cytochrome b's. Hemoglobin and myoglobin 
reversibly complex molecular oxygen. Peroxidase and cata- 
lase oxidize a variety of substrates with hydrogen peroxide. 
The function of the cytochromes involves oxidation and 
reduction in the electron transport chain. When the metal- 


Loporphyrin is a functional group in -a protein it undergoes 
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co-ordination with some amino acid residue of the protein 
and the solvent. In ferrimyoglobin (metmyoglobin) the 
fifth and sixth co-ordination positions are occupied by an 
imidazole nitrogen of a histidine residue and a water 
molecule (Nobbs et al., 1966). In ferrous myoglobin the 


water molecule is absent. 


Properties of Hemin in Solution 

Hemin, when dissolved in an alkaline solution, 
is known to form dimers and polymers (Clark and Perkins, 
BPo4OreCowg1t! and Clark, 1952; Davies; 1952;. Inada and 
Shibata, 1962; Maehly and Akeson, 1958; Jordan and Bednarski, 
1964). When the pH is lowered in aqueous solution the first 
change is protonation of an hydroxide ligand which is 50% 
complete at pH 6.7 (Shack and Clark, 1947). It is postulated 
that this conversion is from the dimeric diaquo species to 
the dimeric hydroxyaquo species. The formation of the 
dimer is known to be inhibited by large concentrations of 
ethanol (Davies, 1940; Maehly and Akeson, 1958; Vestling, 
fo20npecthyhenenglycolm( Sin th folosi) fepynidines(Phiimipse 
1960) and dimethyl formamide (Scheler, 1960). Aqueous 
detergent solutions are also known to favour formation of the 


monomers (Dempsey et al., 1961). 
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Equilibria of Ligand Substitution on Hemin 


A great deal of data on the binding of ligands 
to hemin in a variety of solvents and solvent mixtures has 
been published. Studies of the systems most relevant to 
the imidazole binding study are reviewed here briefly. 
Shack and Clark (1947) have studied the reaction 
of cyanide, pyridine and the acid ionization of the dimeric 
hemin in aqueous solution. Both cyanide and pyridine 
cause the dimer to break up to form the monomeric disub- 


stituted complexes. The equilibrium is of the type: 
Bas 2ML.5 Gls 05)) 


: where M, is the dimer, M the monomer, and L the cyanide 


2 
Swapyri dineseiigandsvsetCowoilideand Clarks (1952)cstudied a 
series of substituted imidazoles in co-ordination with iron 
(III) mesoporphyrin and found a rough relationship between 
the acid dissociation constants of the imidazoles and their 
equilibrium binding constants. Scheler (1960) studied the 
binding of imidazole, cyanate, azide, fluoride, and 

formate to hemin in dimethyl formamide-water solution and 


found that for formate, fluoride, cyanate, and azide the 


éguilibrium jwas of the type: 
M+L — ML 0G) 
and for imidazole and cyanide of the type: 


M + 2L en ML» (07) 
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Gallagher and Elliott (1965) studied the titra- 
tion of hemin with pyridine over a wide range of pyridine 
concentrations and obtained evidence for the break up of 
the dimer to form the monomeric pyridine substituted com- 
pound. A similar effect was noted with caffein and imid- 
azole: (Galtagher®"and ELl-otty 1967; -1968) “in°which “they 
proposed that the dimer was formed by hydrophobic inter- 
action between the two planar aromatic rings of the porphyrins 
leaving all four co-ordination positions of the dimer open 
to the solvent. A spectal study of a histidine anda 
histidine-glutamic acid copolymer complex with hemin in- 
“dicated that at pH 6.2 one co-ordination position of the 
hemin was occupied by the imidazole of a histidine residue 
and the other by a water molecule or a glutamic acid 
residue. On raising the pH to 7.2 a dihystidine type 
spectrum was observed. The binding of a lysine polymer 
with hemin gave a compound that was spectrally analogous 
Weclre rere soycocirome *c"(Brauer “atid Ehrenberg,’ 1963)’ “Ini- 
idazole, pyridine, and glycine binding to the ferriheme | 
peptide Of cytochrome '"c indicated Tlie oni zation of an 
iron co-ordinated imidazole with a pK of 11.1 (Harbury and 


Loach - 41.960) 


Kinetic Studies of Ligand Binding to Hemin 
The binding of imidazole to hemin in an ethanol- 


water solvent is described in detail in Chapter 4. 


dathius-op4 
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Riga tii eu ol. ,) 2907; FAngerman et al., 1969). Gallagher 
and Elliott (1968) followed the break up of the hemin dimer 
to form bisimidazole hemin monomer by means of a stopped 
flow apparatus. Fleisher St di, (1969) studied cyanide 
binding to iron(III) hematoporphyrin in alkaline aqueous 
media. A discussion of these results is included in 


Section 423504. 


1.03 Study of Fast Reactions 
Introduction 

The flow technique for the rapid mixing of two 
reactant solutions was first developed by Hartridge and 
Roughton (1923). In the stopped flow method two reactant 
solutions are mixed rapidly (vl msec.) and which then flow 
down to an observation chamber for spectrophotometric study. 
The flow is then stopped suddenly and the reaction may be 
followed spectrophotometrically for several milliseconds 
after mixing. Both Chance (1963b) and Gibson (1954) were 
instrumental in bringing about technological improvements 
that reduced the mixing time to about 1 msec and greatly 
increased the ease of operation and the economy of mixing 
the reactants. A description of a stopped-flow apparatus 
built in this laboratory ts) Géscribed @nechapter 2.” “The 
relaxation techniques described below, also have greatly 


extended the time range for measuring chemical reactions in 
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solution and these include: temperature-jump, electric 
field jump, pressure-jump, and ultrasonic absorption. 

These techniques were developed and refined by Eigen and are 
reviewed thoroughly by Eigen and de Maeyer (1953). The 
minimum measurable reaction half-times, ty, for the various 
flow and relaxation, as well as magnetic spin resonance 


techniques are compared in Table (1.04). 


Relaxation Methods 

When a reaction mixture is at equilibrium and 
some external parameter that determines the equilibrium 
conditions is changed, the reactants "relax" to the new 
equilibrium with a characteristic relaxation time. The 
external parameter may be temperature, pressure, or 
evectric tera, the welaxation time, t, that characterizes 
the valaxaten uo equilibrium is a.functionof the rate 
constants in the reaction under study. The method can be 
best illustrated by deriving a relaxation time for a simple 
chemical reaction: 


a 


ie SOT e Mewes k® CE S08) 


Ky 


for which the differential rate Law is: 


d{c] _ 
- Ser = ke] - kA) (1.09) 


Consider that upon a stepwise perturbation of some external 


parameter the concentrations change by the small amount 
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Table 1.04 


Comparison of Various Kinetic Methods Giving 


Minimum Reaction Half-times of to 
2 


Method ae 
(sec) 
Flow 
Continuous flow eon 
Stopped-flow 10s” 
Accelerated flow ie 
Relaxation 
Temperature-jump ie 
Pressure-jump a 
Electric field jump Lone 
Ultrasonic absorption TOR] 
Magnetic spin resonance 
N.M.R. 1m 
ETSeR. 1d 


2 from Caldin (1964a) 
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AIA], AIBT,* and. A[C} and 


[A] = [A] + A[A] (17105 

[BJ = *(B)-+ ALB] Cet | 

Sere ee: [OF Ciety25 
where the bars represent equilibrium concentrations. From 


the stoichiometry of the reaction it is also known that: 
AE] = --AfA] ‘= <-A([B] (153) 


pubstiCucaig equation (1.10—l.32) an-equation (1.09) 


GA [C] 








pane = {k, ({A] + [B]) + k_j} A[c] 
—.-°> a 2 
- k, [AJ [8] + k_, [Cc] ~ k, (A[C}) (15.014) 
At equilibrium: 
dic. = 
AG =e Gi 262.5) 
and from equation (1.09) 
Kee bCh ake ADB) =O (1.16) 


Also as the term -k (A[C]) 7 is small it can be dropped from 


equation (1.14) and along with equation (1.16) the result is: 


= Khe (i Alet fe) ) an plesies: 





Integration of equation (1.17) gives: 


yee Red Con Pt i (1.18) 


20 









Th 
Fa h © 7 " . renee Pa “ = i. 7 ag * eo , 7 
ast ero (aL. .) ea ups nitw pole bos (81.1) ‘ tée 


where ATC], is the concentration change at time zero. The 


reciprocal relaxation time is given by: 


cee (Al oe TE) ec eikey 


1 
Thus it can be seen that the reversible second order-first 
order reaction (1.08) is described by a simple first order 
process. Measurement of t at different concentrations of 
CiAl] +806B]) allows7the determination of themtworrate con- 


stants k, arid k 


o] ie 
Temperature-jump Method 
In the temperature-jump method, the temperature 


of a reaction system at equilibrium is increased rapidly, 


and the rate of relaxation to the new equilibrium is 


followed. The change in equilibrium depends upon the thermo- 


dynamic equation: 


dinke =o SrvAn 
ap aie Cee 0) 





Henn etOr tie reaction under study is zero, its equilibrium 
will not be changed by an increase in temperature. In the 
temperature-jump apparatus, built in this laboratory, 

and described by Hasinorfl et al, ) (1969) and Ellis (1968) 
as well as in Chapter 4, a voltage discharge of 20 kV was 
sufficient to raise the temperature of 1 ml of solution 
approximately 6°. The apparatus is described schematically 
in Figure (1.04). The energy used in heating is equal to 


MevoN where C is the capacitance of the high voltage 
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Figure 1.04: Schematic of Temperature-jump Apparatus 
A: Light source 

B: Monochromator 

C: Temperature-jump cell 

Dee Photo multiplier 

E: Oscilloscope 

Eo SpaLrtigap 

G: High Voltage power supply 
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Capacitor and V is the voltage to which the capacitor is 
charged, if losses in energy due to the inductance of the 
apparatus are negligible. It can also be shown that the 
temperature rise in the solution is determined by the 

total resistance and capacitance in the discharge circuit 
and that the relaxation time of the heating time constant 
for the heating of the solution in the cell is given by RC/2 
(Amdur and Hammes, 1966). The heating time constant for 
the apparatus in this laboratory has been calculated to be 
BDOuUte MSc Cm(ULltsS, 19056). [nus reactions occurring in 
times much larger than 5 usec are considered to undergo an 
instantaneous temperature perturbation. The resistance, R, 
and hence the heating time constant is kept low by the 
addition of an inert electrolyte such as KNO. to the 


reaction mixture. 


Study of Fast Reactions in Biological Systems 

Both flow and relaxation techniques have been 
applied successfully to the study of biological systems. 
These techniques are well suited to the study of enzyme 
reactions as these reactions are often extremely rapid. 
The application of flow techniques to biological systems has 
been thoroughly reviewed by Chance (1963a) and Roughton (1963). 
Hammes (1968a) has recently reviewed the application of 
relaxation techniques to biological systems and Eigen and 
Hammes (1963) have reviewed the study of elementary steps 


in enzyme mechanisms. 
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CHAPTER 2 
Design Improvements to the Stopped-Flow Apparatus 


A stopped-flow apparatus was constructed that 
offered a number of improvements over a stopped-flow 


apparatus in previous use in this laboratory (Ellis, 1968). 


2.01 Design Considerations 

The optical efficiency of the stopped-flow 
apparatus was greatly increased by a more efficient light 
collecting lens system. A conical quartz lens (Esco pro- 
MiCreeOptTealent. cjuartz (Leng, 1/2" OD x.3/4" long with a 
focal length of 1") was placed at the entrance port of the 
observation chamber and likewise an identical lens was 
placed at the exit port. Light incident on the entrance 
lens from the monochromator projected an image of the mono- 
chromator slit on the active surface of the photomultiplier. 
With this arrangement very little of the light emanating 
from the monochromator was lost. 

Cavitation occurs in a rapidly flowing fluid when 
there is a local pressure drop across an orifice or sharp 
edge causing the liquid to vaporize (Chance, 1964). Dis- 
solved dir in. the solution may diffuse intoithe cavity in 
the fluid, forming small bubbles that persist and interfere 
with the spectrophotometric observation of the solution. 
Cavitation may be eliminated or minimized by: 


(i) designing a mixing and observation chamber that avoids 
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sharp edges, turns, and has orifices decreasing in cross 
sectional area as the fluid flows downstream. 

(ii) applying a back pressure to the fluid so that the drop 
in pressure never falls below the pressure required for 
Cavitation (Chance, 1964). 

(111) degassing the solution so that bubbles do not form as 
ae -result Oo: Cavitation. 

The mixing and observation chambers in this 
apparatus were designed to eliminate cavitation by (i). The 
performance of the apparatus was such that neither conditions 
(ii) nor (iii) were necessary for successful operation of the 
stopped-flow apparatus. The channels in the mixing block 
were designed such that when going downstream at any turn 
or orifice a decrease in cross sectional area occurs which 
increases the pressure at these points and prevents cavita- 
tion. Also, where possible, the angles at which turns 
occurred were made obtuse. This has the effect of decreas- 
ing the cross sectional area at the turn to a minimum. The 
previous stopped-flow apparatus did not incorporate these 
features and was subject to cavitation. 

If efficient mixing is to occur in a stopped-flow 
apparatus the velocity of the fluids mixing must be above 
a certain critical value such that the flow is turbulent 
ana eddying occurs? Thais criticaw velocity; Yar in cm/sec, 
has been found empirically to relate to the diameter of the 
tube through which the liquid is flowing, d in cm, and the 


density, po in gm/om> and the viscosity, n in poise, of the 
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liquid by the equation: 


Np n 
U. = pda @2.01) 
where Np is the Reynolds' number, a unitless number equal 
eo-about 2000 (Caldints (e645) 7~in the observation: chamber, 


which is a tube 0.2 cm in diameter, it can be calculated 
that Uy has a value of 100 cm/sec. The experimentally 
measured value of the flow velocity in the observation 
chamber was 310 cm/sec as discussed in Section 2.02. MThis 
value is much larger than that calculated for Uy, from equa- 
tion (2.01) so it can be assumed that the mixing in this 
apparatus is aided by the attainment of turbulent flow con- 
ditions. 

The mixing chamber of the stopped-flow apparatus 
Pomeownitmbag 11 DrinciplLe to that constructed by Ellis (1968). 
In effect, mixing is accomplished by a 12 jet mixer. The 
mixing process can best be described by referring to the 
diagram of the pees ecnamveruin Figure (2.01).° Bach’ of the 
channels leading from the driving syringes split the stream 
Of liguid into@lourtseparate streams) by sthe, "splatter". 

Each of the four newly created streams from the splitter 

are recombined at the extremities of the mixer creating 8 of 
the 12 jets bringing the number of streams down to 4. These 
four streams are again joined in the central portion of the 
mixer giving an effective total of 12 jets for the whole | 
mixer assembly. From the central portion of the mixer the 


single recombined stream leads to the observation chamber. 
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Design Details of the Stopped-Flow Apparatus. 


Figure 2.01: 
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Flow apparatuses constructed by other workers have anywhere 
PLO 2° (otk Lt cmatter, 1.964) tole jets (Chance; 1964b) 

The previous model of the stopped-flow apparatus was driven 
by pushing the thumb on a plunger against the driving 
Syringes. It was found that upon pushing the plunger the 
whole apparatus would move slightly out of the light path 
causing an erroneous voltage reading on the oscilloscope. 
Also it was not possible to obtain reproducible driving 
times, and dead times when operating the stopped-flow 
Manually. When studying very fast reactions it is important 
to have reproducible driving times. 

To drive the syringes containing the reactant 
solutions, a pneumatically operated piston device was con- 
SIeiuCtea andethis 1S sshown aoeeroqures (2.02). “The Laboratory 
compressed air line in series with a reducing valve and 
gauge, set at 15 p.s.i., provided adequate force on the 
driving syringes. The mixing block and driving platform 
are joined only by flexible Teflon tubing and hence vibra- 
tions or movements in the driving platform, are,;not itrans— 
ferred costre optical “path of the mixing block: 

The temperature of the apparatus was thermostatical- 
Ly «cOntrolledsby a Rigen water jackets and baffles. 
Thermostatic control of the mixing block ws provided) for 
by the copper water jackets I and J in Figure (2.01) which 
are bolted onto the mixing block. The driving syringes are 
kept at a constant temperature by their snaivr age water 


jackets completely surrounding the syringes and labelled 
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Parts B to D in Figure (2.02). The area directly under the 
3=way valve (3LL3 Hamilton Syringe Co.), part L in Figure 
(2.02) has its temperature thermostatically controlled by 
a large copper water jacket, part H, bolted to the driving 
PravccorieeireadaLuLion, a Lengescopper helical’ coil’ con= 
taining rapidly flowing water thermostatically controls the 
temperature of the space where the Teflon tubing runs from 
the driving platform to the mixing block. Glass wool and 
aluminum foil were placed around the coil to ensure that 
this gap was filled by dead air space only. 

The apparatus was designed so that the cooling 
jackets holding the Hamilton gas tight syringes (2.5 ml, 
1002LL, Hamilton Syringe Co.) were interchangeable. Thus 
syringes of different sizes may be inserted in different 
syringe jackets so that the stopped flow apparatus can 
function for relaxation concentration-jump experiments. 


(Hammes, 1968). 


2.02 Measurement and Calculation of the Instrumental 
Parameters Describing the Performance of the Stopped-Flow 
Apparatus 

A number of instrumental parameters are used to 
describe the performance characteristics of a stopped-flow 
apparatus (Czerlinski, 1966b). Among them are: 
(i) Mixing time -- the time required to mix two solutions 
in a stopped-flow apparatus, which is approximately 1 msec 


if the critical velocity for turbulent flow has been 
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exceeded (Caldin, 1964c). 

(ii) Driving time -- the time elapsed to drive the syringe 
plungers from stop to stop. This was measured from a trig- 
gered oscilloscope trace to be 55 msec. 

(iii) Volume flow rate -- the rate of volume flow in ml/sec 


Bast any Given point. It is calculated. by dividing the 


volume of solution used in each experiment by the drive time. 


volume flow rate = ml ejected/drive time 


eee oer cert ote ml eer (2.02) 


(iv) Flow velocity -- the linear velocity in the direction 
of flow of the fluid at any given point and is in units of 
cm/sec. It is calculated by dividing the volume flow rate 
by the cross sectional area of the channel, In the observa- 


tion chamber the flow velocity is: 


flow velocity = volume flow rate / cross section area 


Coral em 7se0)77 10.0814 cm” =.3.14x10* cem/sec (2.03) 


(v) Completeness of expelling old solution -- the number of 
times that the observation chamber is swept out by fresh 
solution for each experiment. It is calculated by dividing 
the number of ml used per experiment by the volume of the 


observation chamber 


completeness of expelling old solution = 


0.53 ml1/0.0314ml = 17 times (2.04) 


(vi) Dead time -- the time elapsed for solution to flow 


from the mixing chamber to the observation chamber. Because 
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of the dead time of the stopped flow apparatus the observed 
zero time of the reaction differs from the real zero time 

by an amount equal to the dead time. As it is not known 
exactly where the reactants mix completely as they flow 
through the mixing chamber, or as there is no one single 
point at which this happens, the dead time cannot be defined 
With a large degree of certainty. However, it can be cal- 
culated approximately, by dividing the distance the solution 
travels from the point where the reactants first come in con- 
tact in the mixing chamber to a point one half the distance 
along the observation chamber by the flow velocity in the 
Observation chamber: 


dead time = distance / flow velocity = 


Pome scr Shomer scar 88s gard” of dhes¢219 9) 


There is, additionally, an experimental method 


for measuring the dead time of the stopped-flow apparatus 


and this is described next. Ferricyanide is reduced by 
BecOupucwacid and) the reaction is well. defined. )/It.is first 
Order in each of ferricyanide and ascorbic acid) In addition 


the molar absorptivity of ferricyanide is much larger than 
that of ferrocyanide (Birk, 1969) and both ascorbic, acid 

and its oxidation products are colorless compared to the 
absorbance due to ferricyanide. The experimental procedure 
for measuring the dead time involves studying the reaction 
at a constant concentration of ferricyanide but at different 


concentrations of ascorbic acid so that the pseudo first 
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order rate constant for the reaction is changed. Suitable 
reactant solutions were prepared by dissolving 7.4 mg of 
K 


Fe (CN) LnecouRminoL wacersand lL igm.of ascorbic acid 


4 6 
titrated to approximately pH 9.2 in 100 ml of water. Each 

of the Pee rcane. sien ons were placed in the reservoir 
syringes and the reaction was studied at these initial 
concentrations. The pseudo first order rate constant is 
changed by successively diluting the ascorbic acid solution 
in the reservoir syringe by about one half each time for 
three or four experiments. 

As explained previously the point at which the 
flow has stopped does not mark the zero time for the reac- 
tion but differs from the true elapsed time for the reaction 
by the amount of the dead time. But each of these reactions 
studied at different concentrations of ascorbic acid must 
have had the same initial concentration and hence the same 
absorbance. Thus when a semilogarithmic plot of photo- 
multiplier voltage change, Which is proportional to absor- 
bance and hence concentration at absorbances less than 
0.02 absorbance units) vs. time is made the common inter- 
section of all the straight lines must be at a point 
equal to the true start of the reaction or the dead time 
behind the apparent zero time of the reaction. The plots 
are shown in Figure (2.03) and the dead time was determined 
to be 5 msec from their intersection. The dead time cal- 
culated from the flow velocity was 5.2 msec., equation (2.05), 


and is thus in good agreement with that found experimentally. 
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Pigures2s0she Semileogarithimice ploc of voltage vs. time’ for 
the reaction of ferricyanide and ascorbic acid. The dead 
time, determined from the common intersection of the straight 


lines ;) Ss. Ssmsec . 





The addition of the pneumatic drive to the stopped-flow 
apparatus resulted in a two-fold improvement in dead time 
over the previous model constructed in this laboratory 


Carla Sy, BBG &) . 


2.03 Construction 

Plexiglas was used for the construction of the 
mixing block because it is transparent, reasonably inert, 
and is fairly easy to work and machine. The only materials 
that the reactant solutions come in contact with are glass 
(syringes), Teflon (syringes, tubing, valves), stainless 
steel (tubing adapters) and Plexiglas (mixing block). The 
Secrronal mechnod Of CONStruction (parts A to F in Figure 
(eu t wes. Chosen since it Ofrered case Of construction as 
well as a convenient method of cleaning the mixing and 
observation chambers after a period of use. 

To prevent stray light from travelling through 
thesmixing block, One surface of part D of Figure (2..01) 
was painted black. Tests with solutions of very high opti- 
cal absorbance in the light path indicated that the amount 
of stray light was negligible. 

Previous models of temperature-jump cells and 
of a geome flow apperatus (Ellis), 1963) with quartz 
lenses were found to leak solution at the interface of the 
lens and the Plexiglas whether they were glued with epoxy 
cement or a mixture of Plexiglas and solvent. The problem 


was particularly severe after the apparatus had been in use 
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for several months. The problem was remedied satisfactorily 
on the design of the stopped flow apparatus described here 
by employing a press fit of the quartz conical lenses. 

Parts A and F in Figure (2.01) contain a rubber gasket 

that evenly freerioutes pressure on the lens ensuring a 


tight fit onto the ends of the observation tube. 
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CHAPTER 3 


The Kinetics of the Oxidation of Ferrocyanide by 


Horseradish Peroxidase Compounds I and II 


4701 THereduction 


Horseradish peroxidase (EC 1.11.1.7; donor: 


H50, oxidoreductase) was first observed to form a spectro- 


scopically distinct compound, HRP-I1", With HOO} ever. 30 


1 


Dye 


Abbreviations used: HRP: horseradish peroxidase; 
Pee Petite ht Re ie, Re -@CN, HRPSE: compounds: 15° LI 
and III, cyanide and fluoride complexes of HRP; R.Z.: 


Pavey rrulber ss Polis, HP-biy, HoP-ii, HoP-Il: states of 


2 6 


ionization of HRP-II; P-I, HP-I: states of ionization of 
HRP-I; ferrocyanide: hexacyanoferrate(II)- without regard _ 


to state of protonation; ferricyanide: hexacyanoferrate 


3- 


CAKE Fy: Ky Kp? the eciad dissociation constants -of HFe (CN) ¢ 


a Ea 


tion constants either on HRP-I or HRP-II as indicated; 


and H,~Fe(CN) 


2 D1 etc: successive acid dissocia- 


i and k 


Iobhs Robes pseudo first-order rate constants for 


the reactions of HRP-I and HRP-II with ferrocyanide; 


k apparent second order rate constants of 


Zapp’ “3app' 
the reactions of HRP-I and HRP-II with ferrocyanide; ki 
Kor etc: second order rate constants; v: anitwalsreaction 
Ve OCH EV amare sen Jo: molar concentration and initial molar 
concentration without regard fOestate of protonation; Ver 
ar Aur poe voltage and absorbance at time t and time 


TntLAley ys OV; OA: |v, = Via| cand [A, aa 7 ty yo! reac- 


tion half eiitew uy LOnLc strength. 
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years ago (Keilin and Mann, 1937). MTheorell (1942) 
identified another compound, HRP-I, which was formed prior 
to HRP-II when HRP was treated with H,0.- 

The generally accepted reaction sequences and 


stoichiometries for the reaction of HRP and hydrogen per- 


Oxide with a substrate such as ferrocyanide are: 


k 
HRP + H,O, wel OP Dare) RET (3.01) 
Aehe es 72 
HRP—Tor Fe (CN] oe nee tT + Fe (CN) ¢ (et) 
nee as | a= 
HRB = Tiana (CN) pee ee Tt Fe (CN); (3503) 
(George, 1952, 1953a; Chance, 1952c). Thus HRP-I con- 


tains two oxidizing equivalents more than HRP, and HRP-II 
contains one more equivalent than HRP. No particular states 
of ionization of the reactants are implied. The rate con- 
Stants are lPabelled)apparent |rate “constants, ask TT 


2app 
k , at least, show a marked dependence on pH. Ferro- 


3app 

cyanide was chosen as a substrate because it undergoes a 
Simple one-electron oxidation to ferricyanide without the 
production of free radical intermediates. A recent study 
en the HRP catalyzed oxidation of lumrmnol by H,0, ator 
8.0 has shown that the free radicals produced from the 
substrate are oxidized at kinetically comparable rates by 
HRP-I (Cormier and Prichard, 1968).. Even though ferro- 
cyanide is probably not a substrate of physiological sig- 
nificance for HRP, the analysis of the pH dependence of 
the Kapp and Leiaie data can reveal information about 
1onizgations at the active site of HRP-I and HRP-II. 
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Previous reports have indicated both that the action of 
HRP and H,0, on substrates is pH dependent (Getchell and 
Walton, (701; Balls" and*Hale? 1934. "wrider;"'1962) Pand “pH 
independent (Chance, 1952d). A previous kinetic study of 
the binding of fluoride and cyanide ligands to HRP indi- 
cated the presence of ionizable groups at the active site 
erunrorca ana Alpercy, 19673 BLiis’ and’ Dunford; 1968) 2 CA 
complete listing of the pK's of ionizable groups found on 
PRE are l1rsted in Table (1.03). 

This study shows that the kinetics of both the 
reactions of HRP-I and HRP-II with ferrocyanide show con- 
Siderable pH dependence (ko, changes’ by avfactor of 70 


pp 
and SENS DY ea fraccor JOLr over 50,7000" over* the pH Ganges 
the reactions were studied) from which information’ has 
been deduced about kinetically significant ionizations at 


the active site of the two enzyme intermediates. 


3.02 Experimental 
Materials 

Purified, ivyophilized HRP, obtained, £rom 
Boehringer-Mannheim was further purified by elution from 
a Sephadex C-50 (CM) column. The procedure is described 
in Appendix 3. .The R.Z. of the HRP samples used, measured 
by taking the ratio of absorbance at 403 my and 280 mu, 
varied fromeaco tO 2597 An Ri caot about, 3.0 is generally 
taken to indicate a highly purified sample of HRP. The 


concentration of HRP was determined spectrophotometrically 


#F 
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ataS0d. muse Using) aaemolerm absorptivitye ofg9.ile x 1a me tom > 


(Keilin and Hartree, 1951). 

Doubly distilled water was used in the prepara- 
ELOneOLeaaliesolutions.<» Thestonici strength; of) all” meaction 
mixtures was kept constant at 0.11 with 0.01 contributed 


by the buffer and the remainder by KNO, and ferrocyanide. 


2 


The solutions of ferrocyanide were prepared less than one 


hour before use from Baker reagent grade K,Fe(CN) *3H,0 


4 6 


under diffuse lighting and were protected from light dur- 
ino storage and use. This precaution was taken to prevent 


either thermal or photodissociation. The molar absorptivity 


of ferricyanide was._measured to be 1060 M “cm + at 420 mu 


; i =] 
in agreement with the published value of 1050 M Tom 


feer., 1969). The®*molar absorptivity of ferrocyanide at 


420 mu is much less (%5 mM tom *) Chane that Or ferricyanide 


and can be neglected (Birk, 1969). 


The concentrations: of ‘the H,05 solutions were 


analyzed spectrophotometrically by an adaptation of 
methods used by Ovenston and Rees (1950) and Ramette and 
Beandrord (1965). Pive mil of freshly prepared 0.2 M KI 
solueron, 0.1 mi of-0.5%3 anmonium molybdate and up: to 4 mi 
Of Stock HO. solution were placed in a 10 ml volumetric 
flask and made up to the mark. The molybdate catalyzes 


by H-O a wis. ch. “then: zeacts ) to 


Z Zion 


form an equilibrium mixture of I, 


the oxidation of 21 =toel 
A blank solution was 
also made up in the same way as the sample solution but 


omitting the HO. The blank solution is necessary 








because of small amounts of oxidizing agents present in the 
water that interfere with the determination. The absor- 
bance measurements were made against the blank at 353 mu 
using an apparent molar absorptivity of 25,900 Maem = 
for the aE present in solution. Concentrations of stock 
H,0, SOLUETONS VOLE about.5 “x Lone M can be measured quite 


accurately with this procedure. 


Stopped-Flow Kinetic Experiments 

Kinetic measurements were made on a pneumatically 
driven stopped-flow apparatus which is described in detail 
in Chapter 2. The driving syringes and reaction chambers 
were kept constant at 25.0°. The lamp source and power 
supply as well as the spectrophotometric detection system 
have been described elsewhere (Ellis and Duntord,\ 1968). 
The reactions were followed spectrophotometrically at 425 
mu. The total absorbance change was kept small (AA < 0.02) 
so that the relative voltage changes observed on the 
oscilloscope, AV, were proportional to AA. 

The buffer solution was not mixed with HRP prior 
to the commencement of the reaction in the stopped-flow 
apparatus to ensure that the HRP did not decompose at 
either high or low pH. The rate of the reaction under 
study was always much faster than the rate of decomposition 


of HRP at these pH extremes (Maehly, 1953). 


4] 
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Steady State Kinetic Measurements 

Steady state kinetic measurements ee conducted 
spectrophotometrically on a Cary 14 Recording Spectro- 
photometer equipped with a slide wire for the absorbance 
range 0 to 0.1. The cell compartments were kept at a con- 
Stant température of 25.0°. The reaction was followed by 
measuring the increase in absorbance at 420 my due to the 
production of Gee ape an det Hern ocyanide ; HRP, H5051 
and sufficient water equal to a constant total volume of 


0.08 ml were added with Hamilton microliter syringes toa 


cuvette containing 2.00 ml of a buffer - KNO, solution. 

The final addition, which initiated the reaction, was made 
Peymetirting the solution with a glass rod on which micro- 
liter amounts of HO. solution had been deposited. The 
initial velocity, v, was computed at an extrapolated zero 
time by measuring the slope of the absorbance - time trace. 
Blank reaction rate runs were conducted without HRP present 
and in no case was a Significant reaction rate measured 


FOr the uncatalyzed H.O,-ferrocyanide reaction. 


2 


3.03 kesults 

Neither ferrocyanide nor ferricyanide in solution 
with HRP were observed to change the absorption spectrum 
of HRP at least over a period of one hour. This experiment 
indicates that no spectroscopically detectable complex is 
formed between either of eiteee eompounds (and HRP ©itgaiso 


indicates that dissociation of cyanide is not occurring 
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from either ferrocyanide or ferricyanide. The buffer 
solutions were chosen so that their effective pH ranges 
overlapped. From the continuous nature of the rate data 
as a function of pH it was concluded that there were no 
ee ae ee caused by buffers. 

The worl zattoneconstants Ki\ and K woof HFe (CN) 37 


A B 


and HFe (CN) 2 were obtained by interpolation of data 
Poonam tonne Strength study (Jordan jand\Ewing, 1962). 


The value of Ky and Kp at an Zonic stmengqth of, 0.11 and 


4 


oe7.02 is#7247xA0 °° M and 6.0 x toa: M respectively. 


The other ionization constants for H,Fe (CN) ¢ are larger 


chap Ot M. 


mimeeice, Of the HRP-Il — Herrocyanide Reaction 

EG was found experimentally that HRP-II could 
be prepared by adding 1.1 molar equivalents of HO. 
pion O,o7molar equivalents of p-cresol to a 10% M solu= 
tion of HRP. About 70% HRP-II and 30% HRP were formed 


initially as measured from their molar absorptivities 


(sent, 19550) lie. function of the p-cresol LS. to 


reduce preferentially the more reactive HRP-I to HRP 


The HRP-II was converted directly to HRP over a period 


of minutes as shown by the return of .the joriginalyiRe 
spectrum in Figure (3.01). An isosbestic point at 411 mp 
was maintained, indicating the presence of only two absorb- 
ing species, HRP and HRP-II (Chance, 1963a). The original 


HRP spectrum was obtained when the reaction was allowed to go 
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‘ iJ alii 


betier-s aav9, Sallo2, viaoaete. nessa sei 


= ’ vy a 7, = 




















art 
ineyoix3ek ae sb inayagaxa® te OR - : 
mie xviols teva neous eLew sank suioe fe 

Gotan tees wild tape? potantagiy, 4 
‘hy: eter De ther ie no rJonyh 688, 


“va pegues etoshte vo ranene 


BP ot nos a Metoew sy sb att wan 
A 8 1 el 
P a = i , 
£5 An row .. (42¢ aa aoe 
5 + : 
athe . Beesve olaak mae 


ot ae i) bon “M19 su Lave 
a ; % } af ns 4 ee et - ea "Oe “ 


oltasi eit: mlichee 


’ 
i x in 4 = . 

' ih x 2 = 
t i] t ’ eo a | 


w-i.f petbos yd Deteqere 


2 ein Lieve Cre 160m ee 
4 me peratia HE 39 AO. 


4a 
Hit moxi Beugassm as YOLeES 


“ot4one? att 1 (abaes \se:098) 


(4h, evasdsners Pd odd etic yilesr sheapleng sout 





‘tie ag 


get 


oe nes 








TIME (min) 





380 | 400 420 3 440 
my 
Prqgure- 3702s "Plot of “absorbance ve. wavelength at the times 
indicated showing the conversion of HRP-II to HRP. The isosbestic 
point at 411 mi indicated the presence of only two absorbing 


species, HRP and HRP-II. 
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to completion over a long period of time or upon the 
aco Lelomeorererrocyanide. ‘Sutticient quantities of HRP-IT 
were present up to 0.5 hr after preparation to study the 
Oxidation of ferrocyanide. 

The kinetics of the HRP-II - ferrocyanide reaction 
(3.03) at a given pH and in the presence of a large excess 
of ferrocyanide over HRP-II are consistent with the dif- 


ferential rate expression: 


GQiaRP— Li. = {2 ei 
ee K 3 ope FHRP Ti) (3204) 


The disappearance of HRP-II at each pH and for each ferro- 
cyanide concentration was a first order process and this 

was proved by obtaining linear semilogarthmic plots of AV 
vVs.-time. From these plots values of k -were obtained. 


3o0bs 
Values of K 3 6bs obtained ae pH 5.90 are plotted as\a func- 
@20n Of the ferrocyanidery concentration in Figure (3.02). 
From a weighted linear least squares analysis, giving the 


S0170° straight line’, in Figure (3.02) the following rela- 


tionship was proven: 


4- 
Ke bs = ee oe ] (37.05) 
The slope, oe. of the straight line with its standard 
Beeoree (521, OoljmohOe Men eacu ene intercept with 


rie se teanuard.err Or eS. 5:00) oe Lee Worccedt whihich’ is 
vero within thes interval of the 95s confidence limits. 
It then follows that the total differential rate expression 


ise 


-d{HRP-II] _ Am 
a K3 appl Fe (CN) ¢ ):{HRP-IT} (3.06) 
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ri 2 
[Fe(CN)4"] x 104 (M) 


Figure 3.02: Plot of k 


ms 
sqpe MS: [Fe (CN) ¢ I at pH5. 90) for 


the HRP-II-ferrocyanide reaction. 





The total concentration of 


Waevaboue 15x 107° 


in a LGetold excess at ‘the 


larger excess at pH values 


M and 


HRP 


second jorder ratey constants, k 


mically vs. 


3 


pH in Bagure (3.03) 


in the reaction mixture 


greater than 4. 


the [Fe (CN) 67] was at least 


The apparent 


lowest pH values and in a much 


, are plotted semilogarith- 


app 


and listed in Table (3.01). 


The simplest reaction that accounts satisfactorily 








.. 
. 
C 


for the.k data is: 
3app 
H3P-II + 
A 
Ky 
HopP-Il + 
HP-LI + 
j 
K, 
Pim iif + 


( 


a 








He (CN 6 dpe 
Fe (CN) 3” k 2 “- 
HFe (CN) 2 k 3 
SSS 
Fe (CN) 4 re 
6 4 
—_— 
HFe (CN) >> k 
} iE eae 
4- 
Ce 4 
Fe (CN) ¢ k 6 rs 
HFe (CN) 37 k 5 
SEE 
4— 
Fe (CN) ¢ k 9 


(ore 0,4) 


The vertical arrows denote fast non-rate determining proton 


transfer xKeactions. 
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From reaction) (3407) 2) follows thas: 
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eZ 6 10 


pH 


Figure 5.0576) Sem Logarithmic plot of Kapp and Maepo vs. 
pH for the reactions of HRP-I and HRP-II with ferrocyanide. 


The error limits are based on an estimated standard deviation 


of +102. 





Table 3.01 
Apparent Rate Constants for the HRP-II-Ferrocyanide 


“Buffer Key: A, acetic acid-NaOH; GH, gqlycanée-HNO 


BCIO*NGOH a. ae cielS HNO) 


P, KH,PO,-NaOH. 


37 


C, NaHCO 
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Reaction at 25>s0e,and u-=.0.11 
pH eaaco Buf fer” pH eee Buffer“ 
(ie lcee-) (Mi geceus) 
eo oe 10. GH pelo Seen 0. M 
Pee Gr 10° GH Pa aiar M 
Bese, Seat Katig® GH 2/5 alates VA pis M 
es Oo 0 a kkOO GH Ge Smee? lox 61:02 P 
eee 3.3 3 10° GH Sap ee Oe ee ie M 
Brey 1.1 x LO” A Se i ae. P 
Bear atisng soho A eres Fig” P 
eis aS. 1 exael OF A TA Sabian tnsdeix HO." T 
Payee 2.8 x 107 A ho ee eae P 
Weyl 1.4 x 10> A 7 9h ae aa . 
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33 M, maleic 


37 NaOH; GN, glycine-NaOH; 





A number of other reaction schemes in addition to that 
of reaction (3.07) were postulated in an attempt to fit 


the eeiss data. Reaction schemes that had HRP-II species 


with none, one or two ionizations did not explain the 


Kapp data satisfactorily. It was possible to explain the 


Kapp data with the use of only one ionization constant, K 


for the protonated ferrocyanide species even though the 


ie 


ionization constant for the doubly protonated species, Kar 
is known. The two acid dissociation constants for the 
ferrocyanide protonated species are separated by almost 
PwOLOK gUnits. ..in en on as experiments were conducted 
only as. low as pH 2.80, 1.6 pK units away from pK,, the 
relative concentration of the HFe (CN) = species was very 
low and as such its effect on the total rate of the reac- 
tion cannot be determined. 

In a reaction involving two or more ionizations 
such as reaction (3.07) molecular rather than group ioniza- 
tions are necessarily implied (Dixon and Webb, 1964). 

Group ionization constants, of which the molecular ioniza- 
tion constants consist, cannot be determined from the 


analysis of rate data alone. If the pK values of the group 


ionization constants are widely separated, the ionization 
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constants measured indicate the true value of the group 
ionization constant. The example of the Seance ake 

dibasic acid given by Dixon and Webb (1964) illustrates 
clearly the distinction between group and molecular 


HonmLaatiron Gonstants: 


R-CHCOOH — K 
a7 CH ,COOH IN ay 
igre se R-CHCOOH™ 
CH. COOH cy bn,coor” (3.09) 
KN R-CHCOOH a 
CH.,COOH | 2X 


The relationship between the molecular ionization constants, 


' ' 
Ki and Ky ena <he group LOonization Sees PSS say Kor 
Kiyr and as are given by the following expressions: 
t 
Ky = Kiy ‘+ Sg Vor Oh) 
and 
ales RoxKoy (gos) 
2 Pas Koy 
Two pathways of reaction (3.07) such as: 
alee i goon 
nee Fe (CN) ¢ ree eis (3:12) 
and 
ean aie 
Paria HFe (CN) & a Gaile) 


are kinetically indistinguishable and are combined together 


in one term’of the expression relating k Louthe sspeci ric 


3app 


rate constants, ionization constants, and hydrogen ion con- 


centrations. The only unambiguous rate constants obtainable 
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from equation (3.08) are ky and ke. The differential rate 


Law Loo reaction (312) is: 


_ da[HP-IT] 


Z 4- 
ar = kK, [HP-IT] [Fe (CN) ¢ ] (3214) 


Bnd estore roaat vom (3.513), ' 63 sis 


—— a k, (P-I1] [HFe (CN) 27] (3.15) 
but 
K, [HP-IT] 
Pear oe ee (3.16) 
Reet . 
and 


[Fe (CN) 27] [H"] 


a 


[HFe (CN) 27] (3.17) 


momar Ons =(s.bo) and (3.1/)" substituted in equation (3'Pl5) 
give: 


kK 
Sea > [Hp- ai 3.18 
at re [HP II] [Fe(CN) , ] (S718) 





which from inspection with equation (3.14) is of an identical 
tormeseand’hencesthe two reactions (3.12 and 3.13), are indis-= 
tinguishable. 

The Sep data in Figure (3.03) sat pH values 
higher than 8.5 falls approximately on a line of slope -l 
which indicates that the contribution of k, to the overall 
rate is insignificant (Dixon and Webb, 1964). After a pre- 
liminary graphical analysis of the semilogarithmic plot 
of I yepe VS. PH OTobt ain approximate values for the ion— 
ization and rate constants, non-linear least squares com- 


puter analysis (IBM Share Library, 1964) was carried out on 


the Kapp data to obtain the best fit parameters for 
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equation (3.08). It is necessary to'use a non-linear 

least squares analysis as equation (3.08) contains more 

than two parameters’. In a preliminary computer non- 

linear least squares analysis (IBM Share Library, 1964) 

for values of the parameters in equation (3.08) it was found 
that the parameter (k. + kK, /K,) did not contribute sig- 


nificantly to computed values of k Consequently in 


Sapp 


the final analysis of the k data (k, + k 4K) /K,) and k 


3app 8 


were Set equal to zero. The values of the remaining para- 
meters obtained from a final analysis, listed in Table (3.02), 
were used to calculate the lower solid line in Figure (3.03). 


However upper limits for the value of (k, + kK) /K,) and 


kg can be calculated to be 1 x 10° tfese and 25 ae 


respectively by assuming that their terms can contribute 


—— 


up to a maximum of 10% of the value of k . Any\larger 


3app 


Given a model 

Y; = £(X. 4, Xs or eC Xin boy bo, See, bi) 

which predicts the value, Y, of a dependent variable, Y, 
where the model f contains m independent variables Xia 
and k parameters ee and given n observations 


Kees) pina Eee ee ka ae 


the program will compute the least-squares estimates Di. 


That is the program will adjust the ae to minimize 


h 
ao eae ji 

Le 
(IBM Share Library, 1964). 
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Tablews.de2 


Rate and Ionization Constants®* Obtained by 
Non-linear Least Squares Analysis for the 


Reaction of HRP-II with Ferrocyanide 











i (248 2e0eajus no tue Secuc 
kK 
(« + 3 ; 0 (< 162 M see w) 
K 
A 
kK 
(* se : (oP he MM sec 
4 K 
A 
kK 
(« ne! ; @imGktOMEERII0° Me kec 
6 K 
A 
“ 57 
ke O.G4< 25 .MbUseci  } 
ay 
Ky (346 24). 2 heen oO? M 
K, (6 Rut 1. Byte 10noM 
K (297,90:.3) Be 107 °M 


¢ Errors given are the standard 
meter estimates obtained from 
squares analysis. 


errors of the para- 
non-linear least 
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SONtCErIbDULLonN to Soaps presumably would have resulted in 


their detection in the non-linear least squares analysis. 


On kinetic grounds an upper limit of 1 x wiv Motsec.” LOL 


(k., =f kK, /K,) is still significantly larger than the value 
of either of two other parameters, (ky + keKo/Ky) and (ke + 


KK,/K,) , vee Ont ry DuEe POs the wvalue sfor k The 


3app° 


large size of the k, parameter has swamped out the effect 


L 
of the (k, is k 3K, /K,) parameter in determining the value of 


Nyapp: 
It is not possible to exclude any of the rate 


constants from the parameters for exceeding the diffusion 


controlled bimolecular rate limit of 10” oe Secu (Amdur and 
Hammes, 1966b). 
Paveuice Of the HRP-1.— Ferrocyanide Reaction 


The race constant tor the reaction of .HRP and 


H,0, tO Lom uRe—) (reaction, 3.01). 215 about’ 1 ox O. ut 


seGhe (Ghamece elo a 3b; Chance et “als, 1967) at least in 


the neutral pH region. Advantage was taken of this known 


rate constant and the value for Kapp determined above to 
study the kinetics “of reaction (3.02). By using a lange 
excess Of HjO, (1 x ee M) over HRP (5 x 10°! M) and 

6 


with only a small excess of ferrocyanide (1.5 x Us Migs 
reaction (3.01) occurs much faster than reaction (3.02). 
It can be calculated that under these conditions reaction 


(3.01) has a half life of about 0.5 msec whereas reaction 


(3.02) was observed to have a half life of tens of 
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milliseconds depending on the pH. In fact reaction (3.01) 
was, under the above conditions, too fast to be detected 
on the stopped flow apparatus which has a dead time of 5 
msec. Thus reaction (3.01) can be considered complete 
before observation of reaction (3.02) is begun. Reaction 
(3.03), the reaction between HRP-II and ferrocyanide, also 
occurs simultaneously with the formation of HRP-II from 
Meeectlon iw .02):. Moreover, the ratio of aoe £0 Sees 
is always 40 or more depending on the pH. Thus reaction 
(3.02) is also separated in time from reaction (3.03) and 
can be studied without interference from reaction (3.03). 
The reaction of HRP-I with ferrocyanide becomes 
faster at lower pH and could not be studied below pH 
3.67 on the present stopped-flow apparatus. . At pH values 
greater than 6, the reaction is sufficiently slow that it 
could be studied in the presence of a large excess of ferro- 
cyanide. Under these conditions the kinetics of reaction 
(3.02) are consistent with the differential rate expres- 
sion: 


_ @{HRP-I] _ . 
= Oka Gy [HRP-TI (3.19) 


Non-linear least squares analysis of the AV-time data in 
the integrated form of equation (3.19) proved the validity 


of the rate law. Values of k ata. Constant pH of 


20bs 
7886 "are plotted in Figure (3.04) and show that: 


: a8 [Fe (CN) 67] (3.20) 


2o0bs 2app 


Thus the complete differential rate law for reaction (3.02) 
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0 4 5 
[Fe(CN) |x 105(M) 


. A 
Figure. 3°04:  ,Plotiot Ko obs vs. [Fe (CN) ¢ ifaieph) 7.865 f07r ‘the 


HRP-I-ferrocyanide reaction. The straight line is the weighted 
least squares best fit line. The slope, k 

Se ae 2aPpP 
errors (676.5)0.2chex LO M sec . The intercept, (-0.20 + 


, with its standard 
0.08) score is zero within the interval of the 99% confidence 
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Table 3.03 


Apparent Rate Constants for the HRP-I-Ferrocyanide 


pH k 
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is: 


_ G{HRP-I] _ 4- ‘ 
= aeee 5 en Ne ] [HRP-I] C3241) 
in practice -an integrated form of equation (3.21) (Frost 


and Pearson, 1961) was used in the determination of Kapp” 


A non-linear least squares fit to the AV-time curves were 
obtained using three parameters: the rate constant, 


k AV at zero time, and a small term for the placement 


2app’ 
of the base line at inifinite time. A complete derivation 


of the equations used and the use of the non-linear least 
squares computer program for this problem is given in 
Pependix 1. To ensure a valid analysis 10 to 15 points 
were taken from each-trace. At each pH, three or four 
oscilloscope traces were used to determine an average 


Kapp” 


In experiments conducted at a constant pH but 
at two different wavelengths (411 and 426 mu) values of 


Kapp were found to be in agreement within experimental 


error, which indicates that only a single reaction is 


being observed. In other experiments also conducted at 


Cons tent) pH, anewhich [H0,] varied) from ix eo to 


es EO0 « M, values of k agreed within experimental 


2app 


error. The validity of the procedure is verified and 


it can also be concluded that H does not react at any 


eS. 
measurable rate with HRP-II to reform HRP-I. Under con- 


ditions where the [H,0,] is larger glows M) than that 


used in this study (io74 M), H,0, reacts slowly with HRP-II 


to form a third compound, HRP-III (Keilin and Hartree, 
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Polar the eyes results apesplotted’ vs. pH in Figure 


(Se.00) Vandsare listed in Table (303). 


The simplest reaction consistent with the k 











Zapp 
data has a single ionization on HRP-I: 
ei (CN) As ero a 
HP=i 45 
pee 10. 
Ky ee ey) 
3- k 
a C1 (CN) 6 I 
P-I + 
ie 4- k 
Fe (CN) ¢ eZ a 
from which 
sh el k, 4K Ione 
9 jae 4 eae 
St ks et + 
Ky 10 Ky (HY 
i, = G32 3) 
app 
Ky [HT] 
[H ] A 


From non-linear least squares analysis of the Kapp data 


in equation (3.23) the values of the parameters listed in 
Table (3.04) were obtained. The upper solid line in the 


semilogarithmic plot of Sess VS. pH ine rogure- (3.03) was 


computed with these values in equation (3.23) 


Seeeady State Kinetics of the HRPseedtaly zed she. 


2—2 


Ferrocyanide Reaction 

The stoichiometry of the H,0,-ferrocyanide 
reaction was checked at each pH value as both the concen- 
tracion. off tie H,0, added and the ferricyanide produced in 


the reaction mixture were known. Within experimental 


error, two moles of ferricyanide were always produced per 
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Table 3.04 


Rate and Ionization Constants® Obtained by 
Non-linear Least Squares Analysis for the 


Reaction of HRP-I with Ferrocyanide 








Pa (ee £062)" s 10° Netsecr- 
Oak 
k + ij'4 (7m 5 +916) 0x Ton i hae 
10M Sank 
A 
k [oAOMROR es 10 Me secs 
12 
eG 
] i YS AOE Sn ean 


° Errors given are the standard error of the para- 
meter estimates obtained from non-linear least 
squares analysis. 
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mole of H,0, added. This check showed that no detectable 
amount of exogenous or endogenous donor reacted, or less 
than stoichiometric amountsof ferricyanide would have been 
produced. 

The initial velocity expression under steady state 
conditions is derived in Appendix 2. Also included in 
Appendix 2 is a calculation of the attainment of steady 


state conditions which show that the steady state is 99% 


MoEmened 2homsec arter the Start of the reaction. 


fen 
a Bilao Cape oapee to 1p [#50514 
4- 
(Ky app* 2app 5 Kiapp*3app’ '#2°21o ce eee ae 6 J 
(3.24) 


where v is d[Fe(CN) 27], /at. The proportionality of v and 
[HRP], was checked at pH 4.79. The values of [Fe (CN) 271, 
and O71, were kept gconstant at 1.00 x HOM Mend 55 3 
107P M and v was measured as a function of [HRP], . The 
results are plotted in Figure (3.05) and listed in Table 
(3.05). The slope and intercept of the weighted linear 
least squares best fit line, with their 95% confidence 
limits are (42 + 4) sec t ahateh AUeee en 2) OSS ome Myaseo ke. 

The intercept is zero within the 99% confidence interval, 


which indicates that no additional constant term is required 


in equation (32:24).: 


From the results above it was found that app is 
at least 40 times larger than Kapp" With the assumption 
ice a and rearrangement of equation (3.24): 


Zapp 3app 
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Figure, 3.05; 


PlOn ot initial steady state velocity, v, vs 
[HRP] , G@teopitaes 2 and Constant initial concentrations -of 
ferrocyanide and H,0, of 1.00 x oe 
respectively. 


Means SO 
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Table 3.05 


initial Velocity Data for the HRP Catalyzed 
Oxidation Of Ferrocyanidqe at,.25.0°, = Ort 


and pH = 4.79 


Sa [HRP] 9 
(M sec 7) (M) 
25x 10. ° prodxad Oa” 
ino wipe dna’ Seaome 
ie 28x D0aR pro xal0ue 
Sect oald & 2A0.x 10. ® 
Ged xeics! 1 Gbxeltec 
de AD xe O 0d 9.8 x LOT” 
44%, 40) | Ha Beale 
lod, welLOS, 2foesndiee 


4- 
The [HO] 5+ and [Fe (CN) ¢ Jo were 


held constant for each determination 


eS cae ede Gk Sek ee Se 


respectively. 
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= Als Tei, Bil ela (3.25) 


AS vapog et 3app 


PeoneeduatLone(s, 20). t aS predicted that. a plot 


of [HRP] )/v vs.1/[Fe(CN) 6] at a constant [H,0,],) should 


0 


be linear with the slope equal to 1/2k and an intercept 


3app 


equal to 1/2k Data obtained at three values of 


Lapp [#2°219- 
Daeare plotted in Figure (3.06) and listed an Table (3.06). 
Weighted linear least squares analysis gave the rate con- 
Steanes and error limits listed in TabVe (3.07). 

hep cansalso bevshown, from equation (3.25) that 
ao Lot Of [HRP] )/v vs. 1/[H,05] 4 at constant pH and [Fe (CN) 6], 


should be linear with a slope of LDS cases and an intercept 


4- : ; 
of a app oh eN) G ] Data obtained at a single pH of 


0° 
iaoewace plotted inf Figure (3.07) and listed in Table (3.06). 
The solid lines in Figure (3.06) were calculated from a 
weighted linear least square analysis. 

Values of k and k obtained froma 

lapp 3app 

weighted linear least squares analysis of the data are also 
ected in, Tablaewter0 /) . 


Values of k obtained from the steady state 


3app 
kinetics and from stopped flow kinetics agree within the 
error limits shown in Table (3.07). Wihe®enmoom\yin the 


determination of k by the steady state method is 


lapp 
necessarily large because of the small value of the inter- 
cepts in Figure (3.06).and the uncertainty in drawing an 


accurate slope through the data points in Figure (3.07). 
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pH 8.73 





Figure 3.06: Plot of [HRP],/v vs. 1/[Fe(CN)? ] 


values’ shown. * For each plot, pH, [H. 


0 at the pH 

05), and [HRP] ) are brs 
constant and their respective values are: 4.82, 3.55 x 10 M, 
Meio hi0G 4524, 3.79 810s” Mand 9.00,.% 207° m: 8.973. 


Or ioc weand. 4.50.x 100! om: 





Tab le3..0G 


Initial Velocity Data for the HRP Catalyzed 


Oxidation of Ferrocyanide at 25.0° and yp = 0.11 


[HRP], ie u . [HRP] , pH 
v Lre (CN) 1, 9 2%! 0 
(sec) (ie (m+) (M) 
6.9 #110°° RUL0 Mec eros |. OD x Mol?) || ape 
peueatom §265 < lcr) See x lof) 2.0 x to? “age 
‘eee aioe Pig dee ccd lo 20 x Nou 482 
eG 10 1 ar bee te le P20 x Aon © 482 
2.0 x 1072 it eee aa | ome fo” «| a8? 
405 2) 10 17 Sip "ee eo | 200 x Oo”) 4 82 
ery sd ON PilOeibee ter. taf OtOe to” CO 82 
i eUbe SiGe 0°} 265 x TO’) 4a5ax ao’ | § & 73 
BET G1 104" So ee ose 0) ee eo 8.78 
OMCs mee * LOC) Meese mo} 48se@x lo.’ 8.73 
eS xl TOTS mein |. 7295 Sto | 445 x fo 4 83©6 88.73 
1.0 x 10° eld 5 0 | 48S x 10! 8.73 
2.0 x 10° me ee fo 8 ee x fl0n, «| 8 473 
3.2 % Lo Pian SiO Bee xelde 8 473 
255 ed LO; rl elo lO-# S80 x elo. 24 
265 fe Por 2iL04 Moho me WO} 90 x elon’ 24 
7.9 st LO) - 0g) 2.6 x10) “oun Mel Re. ed 
1.1% fo) 10/5 BIG De Gwen OF oa bee we 24 
2 0a 0 of) afore 0 aoe 0 ey eo 
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fet seeps 2 107 286 Se to fb Sotcexela ae meet 
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Biguress 0) Lotsa [HRP] ,/v vs. 1/[H,05], ieee Olt eA oO 
and at constant initial concentrations of ferrocyanide and HRP 
offi 400) Ker0 Rida Wisshe por ls respectively. The solid line 
was calculated from a weighted linear least squares analysis. 


The large errors are caused because the term 1/2 k [H,05], 


yLapp 
is small. 
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This arises because the term 1/2k [re (CN) 17] an 

3app bye 
equation (3.25) is, under typical experimental conditions, 
approximately 100 times larger than the T/P2Ky app 299 !5 term 


effectively swamping out the smaller term. 


3.04 Discussion 

The steady state kinetic study has shown that 
equations (3.01-3.03) represent the minimum reaction 
sequences for the HRP catalyzed oxidation of ferrocyanide 
by H505 ateasrixedapH. “ThesinclLusionyolstast .non-rate 
determining steps would also lead to an equation of the 
same form as equation (3.25). Since the values of Saye 
obtained from the steady state method agree with those 
obtained from the stopped-flow studies, no detectable loss 
of activity of HRP occurs in the steady state reaction 
cycle. 

A number of kinetic studies, including the 
present ones, have indicated the presence of ionizations 
on HRP associated with the active site and are listed in 
Table (3.08). It is generally agreed that the pK values 
on HRP measured to be in the range 10.8 to 11.3 are 
associated with the formation of a hydroxide complex (Ellis 
and Dunford, 1969). The temperature jump study of fluoride 
binding to HRP: over the pH range 4.1 to /.9» (Dunford and 
Alberty, 1967; Ellis and Dunford, 1968) showed the presence 
of two ionizations at the active site with pK values of 


4.3 and 6.1. In substantial agreement was a subsequent 
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Table 3.08 


Tie; pk Values for Ionizations inthe Active Site 
Sreine aeduceo| from KinesalG Seudies at 25.0° and u = 0.11 
The values listed for each of the first and second ionizations 


are believed to represent pK values of the same group, with 
changes in the constants reflecting local environmental changes. 


pK Value 
BLS: Second Thar 
Tonization LTOni zation Tonization Ret. 
Species (60h) COG 
P ordination 
position) 
HRP-II Sew Dae 8.6 TRis. work 
HRP-I - 53 = THasswork 
HRP Jgph = wie 10.8 a 
e 

HRP-F Bee ~ aD 
HRP-CN Ga EOS) a 
a 


(Hii tswane Dunreora, L968) . 
(Dunford and Alberty,; 1967). 


Mean of two independent values which agree within experi- 
mental error. 


pK,values,of 4,and 5,were deduced from small spectral 
changes (Theorell and Paul, 1944). 


© The study could not be conducted above pH 7.9. 
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study of cyanide binding to HRP over the pH range 4.2 to 
11.3 which showed the presence of three ionizations at 

the -dctive site with pk wwalues of 4.1,,.6.4;,. and, 10..)8: (Ellis 
and Dunford, 1968). Both the cyanide and fluoride ligands 
would be expected to carry a negative charge to the active 
Site. In both studies it was found that the pK values of 
the molecular ionizations were increased when the com- 
plexed HRP was compared to free HRP, as might be expected 
from simple electrostatic considerations (Alberty and 
Massey, 1954). 

If the mean pK values of the molecular ionization 
constants on free HRP of 4.2 and 6.2 are related consecutively 
to those on HRP-II, then decreases of about 0.8 and 1.0 pK 
units appear to occur when HRP is converted to HRP-II. 

This would indicate that HRP-II has a more positively charged 
active site than HRP (Alberty and Massey, 1954). 

MheenolLeculareronuzabvonyot pk4ds3+0n-HRP=1 
appears «oto icorrelate;with»the»pK of 5.2 on,;HRP-II,; which 
implies that the effective charge on the active site of 
HRP-Ieis ee ree Eromeathak,olsARe-LiaseaNoj;ionization on 
HRP-I was observed that corresponded to the ionizations with 
pPKrvaluescok 3e4consHRP=LLtors pKa4 s25 OneHRP ,ewhichemay 
indicate that it has disappeared. Alternatively, if the 
corresponding pK on HRP-I were shifted to a smaller value 
it might not have been detectable as kinetic experiments 
could be conducted only as low as pH 3.67 on HRP-I in 


this study. Also an ionization on HRP-I corresponding to 
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the one of pK 8.6 on HRP-II does not appear to be present. 
A number of structures have been proposed for 
HRP=I and HRP=-I1,... For HRP-IT) these include the Fe(IV) and 
the Fe(IV)0O structures proposed by George (1953c). Brill 
and Williams (1961) agreed that the ferryl structure des- 
cribes HRP-II but proposed for HRP-I the existence of two 
components (1) a simple complex with hydroperoxide in the 
Spxcn CO-Ordination position and (2) an oxidized porphyrin 
ring compound. Brill and Sandberg (1968) later put forward 
peopesals*® ror the structure* of" HRP=I including’ (1) iron 
Wnvarquedrivalent state plus a free radical or (2) iron in 
the ferric state=plus a biradical on the porphyrin or 
Peeee rm eeMossrecsal: (1909Mstated tChatwthe results” of 
their Mossbauer study on HRP-I and HRP-II were compatible 
with a Fe(IV) type structure for HRP-II but they also con- 
cluded that there was no difference in the configuration 
of the iron between HRP-I and HRP-II, indicating that the 
second oxidizing equivalent on HRP-I is not stored on the 
BLoies Peisachvetwaly (1963) "and = Brumberg et al. “(1968) 
put forth a series of comprehensive proposals for the 
structures of HRP-I and HRP-II consistent with low tempera- 
ture optical spectra, electron paramagnetic resonance 
spectra and previously reported magnetic susceptibility 
data. "Structures “are written “for HRP-If fn “which two 
oxygen atoms reside in the sixth ligand position of the 


iron atom: 
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[HRP (heme ae oon) tY 
enc h ther toni 2ed form 
[HRP (heme coe CON 


por tiesceruclure Of HRP=l le was suggested that the 7% 
axial ligands are not present or that they provide only 
a weak component to the ligand field consistent with: 


5 Ib AL Tt aL 
[HRP (heme d3 72)! cnt e 


The two oxidizing equivalents are stored not on a ligand 
but either on the protein or the porphyrin represented as 
group xX. 

The results are consistent with a number of these 
structural proposals. The shift in the pK values observed 
in, going Erom HRP to HRP-IT (Table *3.06)- are ‘consistent ‘with 
the Fe(IV) type structure for HRP-II proposed by George 
(1953c) which would have an active site more positively 
Caeeded etiam Proce URE “ine results, of Moss et al. *(1969) 
Poratcacing the CcOnziguratvon about the trom of “HRP-ih and 
HRP-II to be the same is consistent with the pK of 5.3 
on HRP-I and 5.2 on HRP-II remaining nearly constant. Any 
structure proposed for HRP=I- and HRP-II that has “an “active 
site more positively charged than HRP is consistent with 
the shift in pK values found from this study. However, the 
proposals “of -Petsach et al.~ (LoCo) are tie most consis cent 
with these results because of deductions that follow from 
a consideration of their structures. The first is the 


Presence: Of van, 10nization On the ligand of the sixth 
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co-ordination position of HRP-II which may very well 
eerrespong 6G the ionization) of pk 8.6 found from this 
study. The second is that HRP-I has lost its axial ligands 
from the iron and hence the lack of evidence in this kinetic 
Bcucy £Or 2onrzations on HRP-F that could correspond to 
Pie pk Or ©.6 Or HRP=11l ar ta the ones of pK 3,4 on HRP-IT 
end 4.2 OM HRP. This may be consistent with the loss of 
ecitais iiganas ~romsHRP-P proposed by.Peisach ef al..(1968) . 
An alternate explanation foreens, lack of kinetically 
distinguishable ionizations on HRP-I might be that the 
reaction site of each of HRP-I and HRP-II with ferrocyanide 
may be different. Consistent with this conclusion is the 
eiiegese1on by Moss et al. .(1969), that one of the -oxidizing 
equivalents on HRP-I is not stored on the iron. If the 
Oxidizing equivalents are not transferred through the iron 
in each case then the kinetics of the oxidation of ferro- 
cyanide by HRP-I may be reflecting a different set of 
ionizations than the HRP-II - ferrocyanide reaction. The 
fact that all of the ionizations appearing on free HRP can 
be correlated with ionizations on HRP-II could also indicate 
that the reaction site of the cyanide and fluoride ligands 
with HRP is the same as the reaction site of ferrocyanide 
with HRP-II, namely at the iron where it has been proposed 
thesoxidizing Ssquivalent of HRP lig stored (Moss .etlal., 


1969; George, 1953a). 
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The reaction of HRP and H,0, co form HRP~I is pH 
independent at least over the neutral pH region (Chance, 
17 Sh eace ore al 6 owen May noLcate—thatthe 
neutral HO, Species takes part in the rate controlling 


step. Hence HO. would not be expected to be as sensitive 
to ionizations at the active site as would the highly 
charged ferrocyanide species in reaction with HRP-I and 


BRE-it. 
Whe pe OF 3.4 On HRP=It and 4.2 on HRP fall in 


the normal range for a carboxyl ionization. Carboxyl groups 
on horseradish peroxidase include glutamic and aspartic 
acid residues and two propionic acid groups of protoporphyrin 
IX (Shannon Sime Wy OG)» pine phi values «or 5,24 On HRP-I1, 
feo On BRe=L, and 6.2 On HRP fall. most closely in the range 
of an imidazolium group of a histidine residue, of which 
horseradish peroxidase has an average of one gm residue 
pe 00 ‘Gmeoteprocein (Shannon et al., 1966)... Similar cor-— 
relations have been made for lactoperoxidase (Dolman et al., 
L960). 

Brill (1966) in a recent review has compiled 
rate constants obtained mainly by Chance and co-workers for 
the reaction of a number of organic substrates with HRP-I 
and HRP-II most of which were obtained primarily at only a 
single pH value. These results are listed in Table (3.09) 
along with the apparent rate constants for the oxidation of 
ferrocyanide at the same pH by HRP-I and HRP-II. The rate 


GOuSsucants Givensin, Table (3.09) are all apparent rate 
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Table 3.09 


Apparent Rate Constants for the Reaction of HRP-I 


and HRP-II with Substrates® 


Kapp 3 app 

Substrate (Mw *sec™*) (msec *) pH 
Ferrocyanide> ox Loe Ta ae Me 7 
Ferrocyanide? oe 10° 208 10° 4.7 
Guaiacol on Toe ein oe 10> z 
Ascorbic Acid oss ene 4.7 
Aniline site Lax ihe i 
Pp-aminobenzoic acid 9 x 10° fe ps 10° —= 
lumino1° eae 7 x 10° g 

a 


Except as marked the data was obtained primarily by 
Chance and co-workers at 25-30° and was compiled by 


Bilis (CL966)). 


Interpolated from rate constants plotted in Figure (3.03). 


From Cormier and Prichard (1968). 
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constants. A comparison of the apparent rate constants 

is not necessarily informative as in the case of the sub- 
strates other than ferrocyanide specific rate constants 
were not determined, Only from an analysis of the pH-rate 
pLolile fomea Subatrate can teésspecificerate constants, be 


obtained. 


‘t 


>=Bt 
ety Ao cape ysyne} x --Adined eho 


7 © 
“2. GAS Si tadve e tO7 elitozg 











’ 
er es 
renee apd acetal tie) jon ai. : 


a 


sot het dere toto eodsxd: “~ ae 
ss 


rid Lbentsriedeb +n Siew 


~ 


. benta: 


78 


CHAPTER 4 


The Temperature Jump Kinetics of the Binding of 


Imidazole to Ferriprotoporphyrin IX 


The material covered in this chapter represents 
a more detailed description of results previously published 


(Hasinoff er ies 1969; Angerman eval. 1969 J. 


ab *incroduce sor 

Ferriprotoporphyrin IX is the prosthetic group of 
a number of heme proteins including methemoglobin, metmy o- 
globin, horseradish peroxidase, and catalase. The simpler 
name, hemin, is henceforth used for ferriprotoporphyrin IX. 
Four co-ordination positions about the ferric ion of hemin 
are surrounded by the approximately planar dianionic porphyrin 
Brags “rhe eferth and sreth \co-ordination posrtions trans 
to each other above and below the plane of the porphyrin 
Pengeare avarrable Mr*solfut bon “Lor vrurther”’ co-ordination by 
erevservent.§ XA-FayYcrystallograpive Sstudrés*(Hoard*eteal., 
1967) have shown that for the five co-ordinate high-spin 
chloroiron(III) derivative of tetraphenyl porphine that 
the iron atom is displaced approximately 0.38 A out~oFr— 
enand of the four nitrogen atoms. A more recent x-ray 
crystallographic eyay (GOURELTYMan etwas 1969) Mor ule 
bis (imidazole) -a,8,y, 6-tetrapheny lporphinatoiron (III) 


chloride has shown that the iron atom can be considered 


Pike 
eal ee 2 
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Six co-ordinate as the displacement out-of-plane is very 
Siight, .0.009 A. By comparison the distances from the iron 
to the imidazole nitrogens were found to be 1.96 and 1.99 
A. The iron-porphine nitrogen distance was 1.99 a. 

In aqueous alkaline solution hemin is known to 
dimerize (Clark and Perkins, 1940; Cowgill and Clark, 1952; 
Davies, 1940; Inada and Shibata, 1962; Jordan and Bednarski, 
1964; Maehly and Akeson, 1958). It has also been established 
that the monomer-dimer equilibrium is shifted in favour of 
the monomer in solutions containing a high concentration of 
ethanol (Davies, 1940; Jordan and Bednarski, 1964; Maehly 
and Akeson, 1958; Vestling, 1940; Davis and Martin, 1966). 
Maehly and Akeson (1958) found that 50 volume % ethanol 
caused the dimer to dissociate completely to the monomer. 
These observations were confirmed by the polarographic 
study of Davis and Martin (1966) which showed that hemin 
undergoes a reversible one electron reduction in aqueous 
ethanol. Thus, by working in aqueous ethanol solutions, 
complications due to dimer formation may be avoided. The 
kinetics of imidazole binding to hemin were studied in 
solutions that were 44.5 weight % in ethanol. 

Imidazole is a part of the amino acid histidine 
and because of its properties it has been implicated in 
a number of enzyme reactions. The imidazole of histidine 
occuples. the firth Co-ordination position, OL the iron of 
both myoglobin and hemoglobin. Ligand binding to heme 
proteins has been used as a probe to obtain information 


about the active site of proteins (Diven CUval yt L905; 
















: | ae + 
‘(ae oe 
wele-to-sha shen iceih ode 26 oJs0 ihro-00 aie _ 


% ; ok 
betes | ‘aqwaea ya A 000.0 ,anpitia 
* eqepoudtia slosebingveds oF 

, | ss -S 
etinuseu-notk edt AL 

BNO PS t I a 

i Bas tte l9D) ssizomi®’ ; 

; = 
wm ebacnI ;OsefL . serve 

22°, ,meeewtA Das vines 7 een 

ai >~reanonon eft > te 


» 


‘ jaLoe,ai .smnoHnom 6aa 


» 
Th: _goivsd) Soasagae 
[guoy +82°4 \povedtk Bae ‘# 
< 6 


76024 Bas videam 
sy tt’ peat b- eis beauas ae 
arow aneisevoendo sagiiim 

| bre akvelowe youd] 

- pdubiananal & asopisbar 

chow ye aad? .  fomamsem 

wif oh awe notes ques | 

‘ced , L730) 02: oa bid siose Stine: 3 eak+onb 
radia a rMdpisw By bs’ sicu 26fd enolsu los 7 


rs . UJ _ 
BIGLADAME “4 eBay « 
ee : Re 















» 


pesteia bios Onims od% bo Stay 6 2a) 


ras 
- ‘ - * <>», ‘ - 


hg a 
sy Ret ca 


ae ee 





80 


Coredsack etsal., 1966; Duntord and Alberty, 1967; Ellis 
aUCmDUttord ~ 69705). LC was tellethiat the binding of 
imidazole to hemin would provide a useful model for the 


binding of ligands to the heme proteins. 


4.02 Experimental 
Apparatus and Method 

The kinetic measurements were conducted on a 
temperature-jump apparatus constructed in this laboratory 
and based on the concepts put forth by Eigen (1963). The 
electronic detection circuit, monochromator, light source 
and power supply were identical to those described else- 
where (Ellis and Dunford, 1968). The temperature jump 
was obtained by triggering the discharge of a high voltage 
Pamacttorm (Plastic Capacitors Inc., LK 500-104 ZNO a Ou piagle | 
50 kV). The trigger module (TM-11) was obtained from E G & 
G Inc. The voltage for charging the capacitor was supplied 
by a Spellman high voltage power supply model 2040. The 
temperature-jump cell was constructed of Plexiglas and its 
geometry was such that 1 ml of solution was heated. Two 
quartz cylindrically shaped lenses provided for the trans- 
mission Of light through the solutionvin tne cell fhe 


cell was constructed in such a way that it was possible to 


perform experiments on a minimum of 5 ml of solution. The 
cell fitted snugly into a brass thermostatic jacket enclosed 
in a light-tight cannister. The temperature-jump cell was 
thermostatican ly  controlied at@l9.0" “andar voltage: of 20.5 kv 


was discharged through the solution which was sufficient 
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to give a 6.0° rise in temperature. 

A monochromator setting of 413 mu was employed 
which corresponded to the maximum of the difference spectrum 
of the two absorbing species in solution. Changes in 
absorbance caused by the temperature perturbation were 
observed by monitoring the output voltage of a photo- 
multiplier circuit on an oscilloscope with a differential 
high gain d.c. amplifier. The relaxation curve was photo- 
graphed from the oscilloscope screen. The single observed 
relaxation time, tT, was determined by the use of a photo- 
graphic enlarger to match the negative of the oscilloscope 
trace with one of a series of exponential curves specially 
constructed for the purpose. A typical trace is shown in 
Figure (4.01). Further details concerning the construction 


of the temperature-jump apparatus are given by Ellis (1968). 


Materials 

Imidazole and hemin of White Label grade were 
obtained from Eastman Organic Chemicals and were used 
without further purification. The solvent, 44.5 weight &% 
ethanol, was prepared from distilled, pid ere ede ee 
water and reagent grade 95% ethanol. The stock hemin was 


prepared by dissolving hemin in 0.1M NaOH and diluting 


3 2 


until the stock was 1.00 x 10 ~M’in hemin and 1.0 x 10 “M 


in NaOH. This stock solution was stored under refrigera- 
tion and from it freshly diluted solutions, 4.00 x 107°M 


in hemin, were prepared for each day's experiments. 


Kinetic and equilibrium measurements were conducted on 
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Figure 4.01: Oscilloscope trace of voltage vs. time. Con- 
centration of hemin was 4.0 x 10° om and Of Amidazole Lo20x 
1079 The™pH was .7.c- and, LOMIC St eceno ri .. wlan. U4: 
contributed by maleate buffer. The horizontal scale is 


1 msec/cm and the vertical scale 100 mV/cm. Monochromator 


setting was 413 mu. 
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tion. For all experiments the imidazole concentration was 
at least’ 50 times Larger ‘than ‘the’ hemin ‘concentration. 
mts the tocal concentration of tinidazole “in’solutiron was 
very nearly equal to the equilibrium concentration of 
imidazole. 

An operational pH scale was established using 
ee method described by Bates et taln™( 1963) Therr ‘nomen- 
clature is used in which the symbol pH* is used for pH in 
a mixed aqueous solvent. An accurate measure of pH* can 
be obtained with an aqueous glass pH electrode in a mixed 
aqueous solvent system provided the electrode is stored in 
eee Olucron Of Composition, “pH “7 "and Yvonre strength simibar 
to that used in the experiments. The glass electrode 
(Beckman long thin probe combination 38183) was stored and 
standardized with a triethanol ammonium chloride-triethanol- 
Sie Curr er ol pits] .05 eBates*et) als, 1963) Ths 
electrode was used in combination with a Beckman Expanded 
Scale ph meter. All solutions “contained a buffer contribut— 
io 9 | ee GS i macs a ln A strength of 0.11 with the remain- 


der contributed by the inert electrolyte, KNO In the pH* 


3° 
range 4.5 tov6.0 acetate butter “was used; 6.2 CO 7.5 matiearce 
burtcer, “G.3 tows. 3 phosphate, buifier Mandec.8 GeO 10S eboracc 


buffer. 


N.M.R. Line Broadening Experiments 
mele eiletle lL « spectra were run on 5 samples of 


hemin dissolved in 48.8 mole % aqueous ethanol in the 


M hemin solutions within minutes of their prepara- 
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3 3 


Concentratsonmrange yl.25 3) L008 to 5.00 x 10 "Mi -The 
line broadenings were found to vary linearly with the con- 
centration of hemin. Shifts were measured from the 
internal reference, sodium 3-(trimethylsilyl) -l-propanesulfonate 
Eastman Organic Chemicals) which was dissolved at the mini- 
mum concentration for the detection of the methyl resonance 
Signal. The line broadening data were obtained at pH* ll, 
but it was found that the line widths were independent of 
hydroxide ion concentration over the range TCs tor GO. lit. 
The n.m.r. spectra for the temperature study 
Pangerman et ual.,.1969) were run on,a,Varian A-56/60 
spectrometer equipped with a variable temperature probe 
and temperature controller. All full linewidths at half 
height and chemical shifts were reproducible to + 2 Hz. 
A Du Pont 310 Curve Resolver was used to separate the con- 


EePourcitongOtstheawaters and the ethanol to.the hydroxy. pro- 


ton resonance. 


4.03 Results 
Monomer-dimer Equilibrium 
A spectrophotometric study was undertaken to 


determine the monomer-dimer equilibrium constant. For the 


equilibrium 
aa 
2MG== =e) (4.01) 


where M represents monomer and D represents dimer, the 


equilibrium constant is given as: 
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[D] 


K = 
[mM] 


1 (4.02) 


The extinction coefficient of the dimeric species 
and the monomeric species were determined in a1.0 x to M 
solution of hemin at a wavelength of 400 mp and at 5 and 
64.8 mole %@ ethanol. It was assumed that the dimer is com- 
pletely formed at 5 mole % ethanol and that the monomer is 
the only species in 64.8 mole % ethanol. The equilibrium 


constant, K for the dimerization, was calculated from 


3 


Tee 


absorbance measurements on solutions 5.0 x 10 ~ M in hemin 


and 48.8 mole % ethanol, the solvent for the n.m.r. study. 


The values for K, at 25.0° and 34.8° are 25 and 34 ut 


for the solvent composition used in the n.m.r. study. MThere- 
fore there is about 10% of the hemin present as dimer under 

ties COnNdab10OnsS ofthe, nm. study .tiuSimiberity, Ky 
measured to be 140 and 520 aa at 25.0 and 34.8° respectively 


was 


under the conditions used in the temperature-jump study. 

The dimer concentration in the more dilute solutions used 

in the temperature-jump study is calculated to be about 0.13% 
of the total concentration of hemin and is insignificant. 
Furthermore at a constant pH* the hemin solutions were 

found to obey Beer's law for a single absorbing species in 
solution over an appreciable range of hemin concentrations, 
indicating that the hemin species did not change as a func- 


tions of concentration. 


The Acid Dissociation on Hemin 


An acid dissociation constant on hemin was 
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measured by spectrophotometric titration. From a pH* of 8 
Bova SOlLucton o.08M in NaOHeno significant: spectral changes 
were noted in the Soret band at 400 mu of hemin solutions. 
Using the absorbance of an alkaline solution as a reference, 
it was possible by using the method of Coward and Kiser 
(1966) to show that the absorbance changes corresponded 


quite closely to that for a single proton dissociation with 


an acid association constant of PK. 6.63. #For the’ equi lis 
brium, 
A + 
M— + MOH + -H (4.03) 


where M and MOH are hemin and the presumed hydroxy hemin, 


ehewacia dissociateron constant is defined as: 


[MOH] [H” ] 


mie ae 7 aM iw (4.04) 


Phe Eitration reswits, plotted as the fraction of the 
PEOeGnated Species present,a, vS. pH* are shown! in’ Figure 
(4.02). Below a pH* value of 5.8, solutions of hemin 

were observed to precipitate upon standing. Thus it was 
impossible to obtain absorbance readings corresponding to 
complete conversion to the acid form of hemin. The scat- 
ter infthe data points in Figure) (4.02). below: pH*35,6.are 
caused by the unstable hemin solution. No attempt was 
made to measure possible pK's occurring in more acid solu- 
tions and kinetic experiments were not conducted below pH* 


6.00. 
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The Acid Dissociation of the Imidazolium Ion 
The acid dissociation constant of imidazolium 
Ton, 


Tie ee oa eee (4.05) 


where ImH and Im are the protonated and unprotonated forms 


of imidazole, is given by: 


+ 
a JOE #. (4.06) 
[ImH ] 
The acid dissociation Constant 7 kK was 


ImH’ 
measured in 44.5 weight % ethanol at a constant ionic 
Serengun Of 0.11 and at a temperatube Gi 25.0°, by titrating 
a solution of 0.01 M imidazole with concentrated HNO, so 


that the volume changeover the course of the titration was 


negligibly small. The inflection point of the sigmoidal 


eetration curve (ayplot of pH* vs .§volume, of HNO, added) 
miaveacecs the acia dissociation ,constant_of imidazolium 
; -7 

On. , Romy! to bel a3 x alo M. 


The Hemin-Imidazole Equilibrium 

When imidazole is added to a solution of hemin, 
large spectral changes occur as shown in Figure (4.03). 
The spectra exhibit four well-defined isosbestic points 
indicating the presence of only two absorbing species in 
solution which, as discussed below, appear to be hemin and 
the bisimidazole hemin complex. No spectral evidence could 
be obtained for the existence of the singly co-ordinated 


species, monoimidazole hemin. If the state of ionization 


« 


88 


. ’ "1 ni 97 scdw 


. a 
COSBDONI 2G) 


cad) will, 


Ae i soa 
; S230 @ 
' - 


ao ee , 
" wh-fiaw. ay 4») a 


rs 





89 


"uTUSU STOZeEpTWTSTq pexetTdwos AT [nF jo 


WHaIgdeds Se St YW SAAND ‘°TT’O Yzbuexzqs OTUOCT Sli pue 9°7 Seine Sie OS T  pur 


OT X¥ [— 70 Bzen yw pue ‘Gg 4D ‘SdBamND ATOZ STOZeEpWT Jo—SuotjzeaquesuGs pue seen aoe 


€- 62 


SEM UTWOY JO UoTJeaRUSOUCD Aoz nw ut -yzbueTeaem *SA SouequOSde so go, amen pee 


(eM 
OS 9 O09 OSS OOS OST OOV OSE 


JOINVEYOS EV 








7 


ers 


IS 


._*. 
*' 





° 





90 


of the various species is neglected the association of 


hemin and imidazole may be represented as: 


Met 2m —-% Mim CA. OFF) 


2 
The corresponding apparent association constant, apparent 


Pecause 19iS a function Of pH“, is given by thewexpression: 








. [MIm2] (4.08) 


The bracketed quantities with bars indicate equilibrium 
concentrations. The proof that the equilibrium expression 
of equation (4.08) is valid was obtained by two methods. 
The slope of a plot of log (RP) / =A)! indicates the 
number of molecules of imidazole binding to a molecule of 
hemin (Scheler, 1960). The quantities Bers and Aare 
the absorbances of hemin solutions in al cm cuvette with 
no added imidazole (reference), a known amount of imidazole 
added (sample) and a large excess of imidazole added such 
that the bisimidazole complex is completely formed. The 
slope of a typical plot is shown in Figure (4.04) and has 

a value of 1.9 which, within experimental error, confirms 
equations (4.07) and (4.08). Also if equation (4.08) 

is the valid equilibrium expression a plot of AA/ (Tm vs. 
AA should be linear with the absolute value of the slope 
equal to the association constant, Saat (Coward and Kiser, 
1966). The quantity AA is (AL-A,) - This type of linear 
plot is useful when a large excess of ligand is required 


for complete formation of the complex. Linear plots were 
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Figure 4.04: Plot of log [(A,-A,)/(A,-A,) 1] vs. log [Im]. 


The concentration of hemin was 4.0 x Lome Moan Loe “was jo 


The slope of 1.9 indicates that 2 moles of imidazole co-ordinate 


per molecule of hemin. 





obtained over the entire pH* range. A typical example is 
shown in Figure (4.05). MThe R spp data asva function of.’pH* 


ere, listed an) Table (4.01) “end plotted in Figure (4,06). 


Kinetics jot the Hemin-Imidazole Binding Reaction 
The overall reaction represented in equation (4.07) 
may be written as occurring in two steps, again neglecting 


states of lonization of the species involved: 
k 


lapp 
Me in, = MIm 
k 
—-lapp 
(4.09) 
Ko app 
Mim im ==) Mim> 
=—2 app 
The relation between the association constant, Rapp’ and 
the apparent rate constants iss 
k k 
_iapp 2app (4.10) 


“app ~ *_japp*-2app 

In all the temperature-jump experiments only one relaxation, 
T, was observed. The observed values of Tt were dependent 
upon the equilibrium concentrations of imidazole and ranged 
from 0.1 to 2 msec. The single relaxation time and the 
lack of any spectroscopically observable intermediate pro- 
vides evidence that the steady state approximation for the 
intermediate is valid: 


d[MIm] 


qt = 0 (4.11) 


Without the steady state approximation two relaxation times 


would. be expected (Alberty et al., 1963). The reciprocal 
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Table 4.01 


Apparent Rate and Equilibrium Constants for the Binding 
Gf imidazole ‘toPlemin at 25° "and: Tonic Strength 0.11 


i oe Kinetic 


Tl photometric DEES 

x 10° Kapp K _2app app 

al OE rey: = See Pe Seger) 
(ME esac 7) x 10 ee) | (sec -) (Migs 
ees oe 5 O76 2.508 470+ 260 2202340 Tae py Ra 
reo 2722609 If 2 cee 580440 330480 ak n0%.o 
Tite EAS esa AAD 16 Bee 7LO250 6404120 Oa! = ORG 
2240 2.14 2.4 =p 650420 505.0 1.4 +0.4 
1.44 +.14 223 eh] 650+70 6204180 Lee) +0.4 
1344 +.14 2.06 £3 760+70 880+ 220 One wan 
eeo 15 220 nae 480+100 904260 Olea 0 ks 
i254 2. 13 Sars ie! 720250 ee OS 0 G5 2S 
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Oeag.s, 09 ue Boi 550410 510+30 ei ess 
OR4- 2.07 ee cae S 550+40 520290 Le 2 
OF 64 £.06 abe ih, ool 440+80 4004210 el = 350 
OF 26,203 O.. 330, = 03 410+20 650260 0.9 £0.22 
OF 125+ . 01 Orie ere Ob 340430 270440 0. 23 "e207 
0.042+.004 On OU S52 200 290230 26210 0.024 204043 
0.0282 2003 0700692. 001 359*100 20230 0. CO42207 011 





Error limits given are the standard error except on the k 


the spectrophotometric K 
app 


lapp and 


data Whichtlare G@stimated at +10%.< 
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Figure 4.06: 
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96 
relaxation time is given by: 


k k 
Zapp lapp urbe 


k 
-2a k 
i. ees lappa (A242) 
ic 
1 + one [Im] 
-lapp 


The derivation of equation (4.12) is given in full in 
Appendix 4. 

Under conditions where the equilibrium imidazole 
concentration is large, equation (4.12) is reduced to 
fe = ky app [Tm Therefore from plots of 1/t vs. [Im] 


fhe Limitang slope: at lavge [am] yields values of aoe: 


Petyprcall limiting slope pfot is shown in Figure (4.07). 


Values of k 


1 obtained in this way are plotted in 

app , 

Figure (4.08) and listed in Table (4.01). Values of 
on and etna’ Se es were computed from equation (4.12) 


using Sad from the spectrophotometric titration, T, 
and [Im] as’ input data for a non-linear least squares pro- 
gram (IBM Share Library, 1964). A complete listing of the 


program as well as its use is given-by Ellis#(1968). The 


values found for the parameters K _dapp and res Sheers 
are listed in Table (4.01). In addition, the three para-— 


MebLers;, Spee! See and a oaee ee 


unknown parameters and best fit values were obtained using 


, were regarded as 


the same computer program. The values of Kapp obtained in 
this way are listed in Table (4,01) also. Additional: proot 
of the validity of the steady state assumption and equation 


(4.12) is the observed behavior of the reciprocal relaxation 
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Figure 4.07: Plot of l/t) vse (im) tor imidazole, binding eto 
hemin at pH* 8.82. The solid line is an estimated limiting 
slope and the dashed line was computed from eq. (4.12) using 


K app determined spectrophotometrically and best-fit values for 


and k 7K 


K_2app Zapp’ “-lapp’ 





klapp x 107°(M'!s-1) 
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Pegqute 4 7UG: PIO OL Kapp vs. pH*. 


that predicted using the best four-parameter fit to equation 


(4720) ; 


The: SOLLG liner as 


The circles are experimental points and error limits 


are estimated. 
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taney, wel /t, aseaatunctionsofs([Im] fasmin Figure, (4.07)).and 
the successful fitting of the 1/t and [Im] data to equation 


F4512) 


The N.MsR. Results 

It was observed that the solvent n.m.r. resonances 
had their natural linewidths when a 100 fold excess of 
imidazole was added to a hemin solution. This result con- 
firms that imidazole binds to the paramagnetic center and 
displaces solvent molecules from the inner co-ordination 
sphere of the metal ion. 

From an analysis of the transverse relaxation 
ee! 55: of the methyl (-CH) and methylene (-CH,-) 
protons of the ethanol as a function of the water to ethanol 
concentration ratio it was concluded by Angerman et al. 
(1969) that the replacement of water or hydroxide by ethanol 


or ethoxide in the first co-ordination sphere of the iron 


may be represented by: 


K 
4 
wen AS 
HOMOH , oF CH, OH HOM(C.H.OH) + H50 (Ayo 3) 


The equilibrium constant for equation (4.13) is given as: 


ahaa [HO] [HOM(C.,H,OH) ] 


4 [C,H,OH] [HOMOH, ] (4.14) 


and was found to have «a value .or G2 eat 35 Si Anger menvet 


aby, L9GoPerRre it is assumed that this vatue of Ky is 


independent of temperature then it can be estimated 


that under the conditions of the n.m.r. temperature study 
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the hemin is 90% in the ethanol form. Similarly under 
the conditions of the temperature-jump study the hemin is 
70% in its ethanol form. 

The analysis of the temperature dependence of 
the line broadening data indicated that the average rate 


constant forpethanol exchangs,is 1.8 x 10° sec and 


for water exchange 6 x for aoc both at 25° (Angerman 
Carel. , 760). fe Wae eiso concluded that 1f a “limiting 


Sy mechanism is assumed for the solvent exchange then 


the mechanism can be represented by: 


ke C,H, OH k_¢ 
mrt ai 5 
HOMOH, <== {HOM! — HOM(G Hon) (4.15) 
k HO k 2 5 
ee) 2 6 
mere lSo totehows thate 
kas 
K, = me = ag ale nic 
5) 6 
: 6 -1 
The values of Ke and Ke are ai 25°) 16 x 10° (sec and 


Om 10° Suess therefore k_¢/kK_s is approximately 25. 


The high ratio of k_¢/kK_s is not consistent with the known 


nucleophilic order of water and ethanol may be due to 
preferential’ solvation of the hemin by ethanol in) the 
outer co-ordination spheres. 


A 04 DIiscussnen 


There are two kinetically indistinguishable 


reactions which explain successfully the pH* dependence of 
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ie an ae . ; 


the ean rate data. In one imidazole and imidazolium 


ion are postulated to react only to the basic form of 


hemin: M 
re Kal ky 
TImH + MOH eee 
steady state 
al intermediates 
Im a MOH rate 
are? 
Ke | ( ) 
M 


The vertical arrows indicate fast non-rate determining 
proton transfer reactions in buffered solutions. The 


horizontal arrows indicate the slower rate determining 


association reactions. According to equation (4.17) the 


apparent forward rate constant is given by: 





= ImH bmnd £4.20) th CLAS) 





ImH* 
KT me || iat 
ao se SS 
is a steady state 
Kall intermediates 
k 
MOH + Im sate} (4.19) 
“i 
ImH” 
for wien Xi app aS 
ie (H’] 
i rgd 


k (4.20) 


K 4 
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Corres FA te a 
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It is easily seen from a comparison of equations (4.18) 
and (4.20) that the equations are of an identical form; 


the parameter kj /K in equation (4.18) corresponding to 


ImH 
the parameter k,/K, in equation (4.20) and the parameter 
ky in equation (4.18) corresponding to the parameter ky 
in equation (4.20). The unknown parameters in equations 
(4.18) and (4.20) were determined using the non-linear 


least squares program that was modified to minimize the 


sums of the squares of the relative rather than absolute 


residuals. Several different approaches were taken in the 
determination of the unknown parameters. In the first, all 
four parameters of equation (4.18); ky, Ky, Kom’ and K, 


were regarded as unknown and best fit values were obtained 
by non-linear least squares analysis. It is readily seen 
from equations (4.18) and (4.20) that the two equilibrium 


Constants, K and bse are distinguishable due to the 


TmH ’ 
symmetry of these two equations. The results of the four- 
parameter analysis are listed in Table (4.02). In the next 


approach, all parameters excepting K for which the 


7 


ImH’ 
valueel,8 x LO ¢.M obtained by titration twasmuseo savore 
regarded as unknown. No ambiguity exists for the deter- 
mination of Ke as existed with the four-parameter fit. 

The results obtained are Listed in Tabke (4.02). 9 The 
discrepancy in the value of Ki measured spectrophotometrically 
and from the analysis of the kinetic data may be due to the 


ethanol and water forms of hemin having different acid 


dissociation constants. The equilibrium could conceivably 
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Non-linear Least Squares Analysis of pH* Dependence 


Rate and 
Equilibrium 
Constants 


for the parameter estimated by the non-linear least 


squares computer program. 


interchangeable giving an equivalent fit to the k 


data as explained in the text. 


Values of K. and K 


Table (4.02) 


or the vk 


lapp 


a, 


Data 


Four Parameter 
Computer Fit 


OA ERP) 
(2 te) 
C2 AB a2 } 
C265725 O21) 
(4.8+0. 8) 


(4.12085) 


ImH 


Rate and Equilibrium Constants Determined from 


AB 4 gS 


Three Parameter 
Computer Fit 


ClO 02) 


AOE OG 


(4,440.2) 


| C2HISA0 22%) 


(or. 8} 


Errors on the parameters are the standard errors obtained 


in the four-parameter analysis are 
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be “of the’ types 


! 
+ Na 
HOMOH ,, + H ——— H,OMOH, 
|| aut (| (4.21) 
ll ina Sa 
HOM(C.H,OH) SH a HOM GOH OH) 


The value of the measured K, would be a combination of on 
and KL’: It is known from the n.m.r. data that water and 
ethanol undergo different rates of exchange from that of 
the hydroxy hemin species (Angerman et al., 1969). 

in edditLion, any reaction that’ is a combination 
of reaction (4.17) and (4.19) will provide an equally 
vValud fit to the *1 app data. 

From an examination of the equilibrium data 


at least two sets of equilibria provide a rough explanation 


of the K data: 
app 


+ ee 
MOBY --) lime ae ljag) —— Mim. + H50 (Ae22) 
£0r which, 
PS a TH) / Ke en Te ee ea a 
app a ImH ImH 
(4. 23) 
and 
iv 
M + 2Im a Mim. (4.24) 
Odw Wik Cry 
pe Ee ait ie ae tries amr pa” (4.25) 
app ImH a : 


From a simple electrostatic consideration the 


presence of an OH ligand trans to the ethanol and water 


104 

















4 : q 
‘ ” : — a 
=. B n 


a — + ms 
OO. TS 2 H+ pie sal ia 


~~. 
1 
i" 
7: 
= 


bi uct ait betvesem af3 to sulsy adg 
- % 


> 
w 


wit mos) awond at #2 a brite ~ 


~ 


7461 20ewitib opreba0 Domenie 


piensa) esicsue oined .vxo thd sdt 


JolLteAer YS ,mor13Lbbs ah: . ~ 


Fu = 
na - 
ee 


Ol.s) Soe (TL.6) sotssss= 


sg Seb .,4 sf3 ef 222 Be 
~ ) qgel ‘- 
.JaatmMexe os mow 
*: / . 7 a 
iilrupe to eyee. ows Geass Aas 


- 
me 
~ 2 ne 


—* sit Og 


6 & 
— aq 


i) 


, a 


fo te MT “wat + FO 
Sue 
. doa tQ2 
ae - | 7 d A 
[ 7 a ‘ ‘7’ y’, + ” ; ‘ 7 
al + 2)” Gat" H} vi, a aga™ 
(28.3) mt Un ee ’ pats a 





105 


on the two hemin species would be expected to have a 
labilizing influence on these ligands. For example the 


water exchange on Fe(OH 50H" is faster than on Fe (OH,) 2" 


5) 
(Ergen andoWwilkins;) e©r965):) 2con this ibasist théevreaction 
scheme (4.19) appears less plausible if the solvent exchange 
rate has an important influence on the rate of substitution. 

Pica Sxl limiting mechanism or dissociation 
mechanism is considered (Seewald and Sutin, 1963) then the 
£rrst reactron vin (4.b7) ibacomes: 


M(C.H.OH) 27 ce Mt?” 4 C_H_OH 
2 aa D725 
b 





(4.26) 


Meee yt es “MIme OH 


where Mi*"is a pentaco-ordinate intermediate and M(C,H,OH) * 
is hydroxyethanol hemin with only the most labile and pre- 
valent ligand shown. The total charges on reactions (4.26) 
are shown as superscripts. Since no inhibition of the 


rate is observed at a large concentration of ligand (Sutin, 


1966), 
ee 
ie Kk, [EtOH] (4,27) 
The value of ie Lop aa 10° eeces From the, Nem.r. resizlts 
and an average k, = 6 x ;iGpeyasie ne eile m+ sect. The 


Correction of 0.7 25 applied to ky since the n.m.r. results 
iicvecatea that. tie hemin 16 /02 "in 1tCs-ethano#® form. The 


concentration of ethanol is 8.2 M so that k 7k, is predicted 
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to be approximately 40. This means that imidazolium ion 

is predicted to be 40 times better as a nucleophile than 
ethanol. Also as ky oe Ky the imidazolium ion is pre- 
dicted to be a better nucleophile than imidazole. Both 

of these predictions are highly questionable. Similarily, 
if the same argument is applied to reaction (4.19) imidazole 
is again predicted to be approximately 20 times as good a 
nucleophile as ethanol - again a questionable prediction. 


However if an ion pair mechanism is considered. 


the faster reaction in reaction (4.17) becomes 


vee ae _ 
j eta ay 
M(CH.OH)° + ImH M(CH.OH) ,ImH 
(2.28) 
g waloes*s qual ts 
M(CjH,OH)* ,ImH” —S> MIm*” + H* + C,H.OH 


y, 

where K is the ion pair formation constant. The ligand 
Substitution) Tateyconst ant isos itiethe,ionj pair equidibraum 
is established rapidly: 


loo 
e 


eee ctr) Sy 72 


As there was no inhibition Of the reaceionvac 


high imidazole concentrations equation (4.29) simplifies £0: 


k, = kX (4.30) 


Also .f£ 20 Ms assumea that the rate of solvent 


exchange is not affected by the ion pair formation, (Sutin, 


1966) then the exchange rate for ethanol is ky = eX 10° 


ite and since kK) = 9x 10° vom eae as before then 
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Beet. Monies The value, of 4 vie is reasonable for an ion 
pair formed by a double and single oppositely charged 
species, as are the ImH* and hemin species. By comparison, 
the ion pair constant between the Co(NH,) .c1** and GHeCOne 
species has been measured directly to be 5 uo (Archer.et 
aloe L965); 

If solvent.removal,is. faster in the ion pair, 
then the value found of K represents an upper limit. On 
@eraind lar basiteeche ion pain’ constant, Lor,.lmyin«reaction 
Bett 1 Se above -()..02 God A decrease in K is expected 
when a cation is replaced by a neutral ligand in the ion 
pair. Similarly, ib can. Lom pai .conoetant is calculated for 


-1 
Boeemuan-a( dal) prabiesuot Kh equakrte 2cand “0.02, M are 


obtained. The K of 2 mt isylLanger,than..that -expected.,foxr 
an ion pairing between a neutral and charged species. Thus 
reaction (4.17) appears the most reasonable on the basis 
Geocnie “Lon pair™ calculation. 

Bleveches et al. (1963) have studied ligand 
binding to cobalt(III) hematoporphyrin IX presumably in 
monomeric form in aqueous solutions. Hematoporphyrin IX 
has two CH.CH (OH) - substituents in place of the vinyl sub- 
stitutuents in ferriprotoporphyrin IX. It can be shown 
that ff. a reaction scheme similarjtc reaction )(4.09) is 
assumed for their cyanide-iron(III) hematoporphyrin IX 
reaction then their observed pseudo first order rate 


constant, ke , corresponds to iar Similar behavior of 


bs. 
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Coyarand Ko are observed at high ligand concentrations 


bs 


for both systems; a plot of Ko Orel /e Ss linear with an 


bs 
intercept of zero. However Fleischer Ou, al .s(1.g00) Were 


unable to observe any leveling off in Ko as they operated 


bs 
at a much larger excess of ligand than in the temperature- 
jump experiments. They were able to obtain evidence for a 
Pm ting Sit ligand substitution mechanism for the cobalt 
HUE noe the TFOnR POrpnyrin. 

The binding of imidazole to ferrimyglobin, a 
heme protein containing the prosthetic group ferriproto- 
porphyrin IX, has been studied by Diven et al. (1965) and 
Goldsack Ctlal.ec.0G)... The apparent essociation rate 
constants were much smaller (< 350 eS than those 


measured in this study presumably because of the hindered 


nature of the environment of the prosthetic group in ferri- 


myoglobin. 
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Appendix 1 


The Analysis of Second Order Kinetic Data Obtained 


on a Stopped-Flow Apparatus 


In the study of a second order reaction such as: 


Ror Bees (Al.01) 


Bor which the rate. law is’: 


=e TA [IB] (Al .02) 


where x is the concentration of B that has reacted at 
time t, it is experimentally advantageous to study the 
reaction if possible by pseudo first order kinetics. 
Under conditions where one of the reactants is in much 
greater excess than the other and remains effectively 
constant throughout. the course of the reaction, it is 
not necessary to know the concentration of the other re- 
actant but only, some physically measurable quantity that 
is directly proportional to it. When the reaction rate 
becomes fast, however, (reaction half time less than 10 
msec) it becomes necessary to lower the concentration, of 
the reactant that is in excess, to slow the reaction. The 
treatment of the data is different since pseudo first 
order conditions no longer exist. 

A method was developed with a consideration for 
the experimental limitations imposed by the stopped-flow 


apparatus that used the integrated second order rate 
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expression (Frost and Pearson, 1961): 


1 In b(a-x) 
(a-b) a(b-x) 








= kt (A1l.03) 


where a and b are the initial concentrations of A and B. 
For the following derivation the following 


assumptions are made: 


Le aaD 

2. A is assumed not to absorb appreciably com- 
pared to B 

3. Voltage changes as specified are proportional 


to absorbance changes which is true if absorbance changes 
are less than 0.02 absorbance units. 

Figure (Al.01) is a hypothetical stopped-flow 
oscilloscope trace where the voltage output of a photo- 
multiplier circuit is measured as a function of time. The 
axes represented by the solid straight lines do not represent 
the true zero time or the true base line at infinite time. 

Starting from the left of the oscilloscope trace, 
the nearly vertical solid line marks the start of the flow. 
The long horizontal portion of the trace is the “driving 
time"; the time period over which the flow exists and has 
been experimentally measured to be about 55 msec. The 
break in the horizontal portion of the trace, where the 
vertical axis is placed, marks the point at which the flow 
is stopped. However it does not represent the true start 


of the reaction as the reactants are mixed not in the 
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Figure Al.01: A hypothetical stopped-flow oscilloscope 
trace of voltagesvss time. Ie solid axes @represent: the 
arbitrarily places axes and the dashed axes indicate the 


true axes at the true zero time and the true base line. 
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observation chamber but in the mixing chamber. The 

time that it takes for the reactants to flow from the 
mMExihg COQ tne observation, Chamber is the “dead time", oe 
and was measured to be 5 msec. These and other parameters 
that characterize the stopped-flow apparatus are discussed 
in greater detail in Chapter 2. The voltage change, AV. 
at the arbitrary zero time does not represent the voltage 
change due to 100% B present. The voltage change corres- 
ponding to 100% B present, AVon is therefore hidden due to 
the finite dead time of the apparatus. 

The parameter 6V represents the error present in 
choosing the true baseline which is theoretically attained 
Beeineineee Emesorenot at all if the product formed decom- 
poses. The horizontal solid line is the arbitrary base 
line from which the quantity AV, is measured and differs from 
the true base line by the amount 6V, which can be positive 
or negative. To each of the quantities AV, and Ae is thus 
addedethe correction iVeed trfollows then «that cthe tiraction 
that the voltage change has been decreased is the fraction 


of x that has formed, hence: 


AV eo 
bas eee se) (A1.04) 


Solving for x in equation (AL 04)« 
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Combining equations (Al.04) and (Al.05) and solving for 


AV, : 


ven uate so PRE waaay. Sy) 
O (9) 


AV (A1l.06) 
Ve (ae (a-b) kt 27h) 


wre equation 1s solved for AV, as this is thé quantity that 
is experimentally measured as a function of time. Equation 
(Fw 06)) ts in a form that is amenable to non-=linear least 


Squares analysis to obtain the three unknown parameters, 


5 AVG. and 6V. A non-linear least squares computer program 
has been described by the IBM SHARE Library program SDA 3094 


(1964) and modified by Ellis (1968) that is suitable to find 


the unknown parameters. 

With equation (Al.06) as the model equation in 
the program and AV, and t asyinput jdata,. initial guesses 
are estimated for each of the unknown parameters k, AVG 
and 6V. With these initial guesses observed and calculated 
values of AV. are compared and corrections calculated until 
the least squares criterion is satisfied. This involves 
minimizing the sums of the squares of the residuals, that 


is the observed minus the calculated values of AV, . 
FORTRAN Listing of Sub-programs and Comments 


The FORTRAN listing of the two subroutines FCODE 
and SUBZ are given on the next page. 
The various FORTRAN parameters and variables are 


defined accordingly: 
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Unknown parameters 


Bel) ao 
Bt 2h Po V 
Bit) ae AV, 


Constant parameters and conversion factors 


Bo) a 
etay: = b 
B(6): the number of graticule divisions per inch of graph 


paper on which the enlarged oscilloscope trace is read out 
on. 
B(7): the number of seconds per graticule division or the 


Sweep speed of the oscilloscope time base. 


Variables 

Y(I): the dependent variable is the number of inches mea- 
sured on the vertical scale of the graph paper and is 
Greccliy proportional to AV, - 

X(I,1): the independent variable is the number of inches 
measured on the horizontal scale of the graph paper and 

Poe pLoportional EO timest. 

The purpose of subroutine SUBZ is to convert the 
experimental measurements of distance on the graph paper to 
true time. This is accomplished with a DO loop which in 
turn, multiplies Gach of the anput data points, x(1,)5, by 
the two conversion factors B(6) and B(7) and to this adds 
the dead time 0.005 sec to obtain the true time measured 


from the true etart of the reaction: 
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Subroutine FCODE contains the FORTRAN version of 
equation (Al1.06), the model equation. Also included in this 
Subroutine ts Che option for printing out auxiliary data: 
the original input data, the true time, and the different 
concentrations. 

In the limit for a >> b, the second order integrated 
rate expression gives the pseudo first order expression and 
the subroutines as written, are also applicable to a pseudo 


first order kinetic analyses. 
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Appendix 2 


Derivation of the Initial Velocity Expression for’ the 
Steady State HRP Catalyzed Oxidation of Ferrocyanide 


by Hydrogen Peroxide 


The reaction sequences of the HRP catalyzed 


Oxidation of ferrocyanide by hydrogen peroxide are: 


k 
cee ea _1apP . ypp-t 





k 
nee Fe (CN) AGN eee Fe (CN) 37 


k 
HRP ar ht Fe (CN) > SIE Serge es Fe (CN) 2 


CAZe OE) 


where no particular state of ionization is implied for 
any of the reactants. From inspection of equations (A2.01) 
the following conservation relations apply throughout the 


eourse of the reaction: 


[HRP] +2 fRRP-T) +. [HRP-leqp= [HRP] 4 (AZ02,) 


[Fe (CN) 27] i [Fe (CN) 27] : [Fe (CN) 67] (A2.03) 


0 


where the concentrations with the subscript zeros are 
initial concentrations. 

For-reactions- (A2.01)* the following differential 
rate expressions apply: 


e315 aay ements 


ye 
at ee tO Boe ig [HRP-I] [Fe (CN) ¢ ] 


2app 


(A2.04) 
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ad{HRP-Ii] _ ee RP-TIilF on) 274 

epee Sosne PS tN) en K3app 4 DEE SEN 
(A505) 

d [HRP] 


Ain 
—3zr CO Reena rs rec oN le ios ee aoc) 
(A2.06) 


The steady state assumption applied to each form 


of HRP is: 


Gipheo= Lise Ure te » GlHRP les 
AES genre O°. Sy lemragsepaereracrns am eagen pase ig 0 (A2.07) 


and when applied to equations (A2.04-A2.06) gives: 


Sy eyo = ee a ] (A2.08) 
Axe eee . Ais 

poate, aims) [eetCN) | ] = UNAS as IT] [Fe (CN) ¢ 1 A209) 
carte ss 

beepe te ttt lee(CNn) | ] = AT app P50) (A2.10) 


Substitution of equations (A2.09-A2.10) in equation (A2.02) 


gives: 
4- 
[HRP-IT] k [Fe (CN) ] [HRP-ITI] 
[HRP-II] + eas + 2 a 
2app Lapp, {i707} 
ae 
[HRP] | (AZ) 


The differential rate expression for the production of 


ferricyanide is: 


d [Fe (CN) 27] a— 
SSS ete [HRP-I] [Fe (CN) ¢ ee We 


4— 
aE Baas [HRP-II] [Fe(CN) 5] 


3app 
(AZ k2) 











a ) oe pT) as , x — St St 5 = — 


14 

a 
< 

ae 

ma 
outs 
- 


an, rei Wh 36 


rr dernw = 
} iz HID , 


“4 ) ob 


so of baiiqas fadw base 


5 @ 
a 


x 7" “ OF ; 73 | Ke 
ca mes ae ak \ eacitel “ws 


» 
(_- 


qos” 


a = (°27(49)89) [lai 


qqpe™ 


’ 


(wo) e%) [11+dnH] 


a - 


solteups t0 aoljusiseceeey 


:eavio 


128 
Substitution of equation (A2.09) in equation (A2.12) gives: 


d[Fe (CN) 37] i 
<t ae tapi 2K 3 app LHRP-II] [Fe (CN) ¢ ] (A2 Uo) 


Solving equation (A2.11) for [HRP-II] and substituting the 


result in equation (A2.13) gives: 


d [Fe (CN) 27] 2k k k [HRP] [Fe (CN) 6] [H,0 


lappyeepp -sapp 2! 
at 


(Ki app 2app a 1 app*3app’ [H,0,] + Seren eae [Fe (CN) ¢ 


(A2.14) 


Defining the “initial "velocity, Py? tas a [Pe (CN) 271, at zero 
time then: dt 
O)4 


A 
+ STapp Capp) 22.0 + eo sae a 6 ] 


ba 
i oy oral ta pars west ds phd SC Dy Jo lH, 


(k 


lapp*2app 0 


(A2.15) 


Both. k and k have been measured as shown 
2app 3app 


in Chapter 3 and it was found that Kapp was 40 or more 


times k larger than k With the approximation (k 


3 app” Zapp 


>> Kapp? equation (A2.15) simplifies to: 
[HRP], 


Vv 


cK he) en a 19° yl wec.” (A2.16) 


yes 
Ge yg HEMT J a app 2°20 


Calculation of the Attainment of the Steady State Conditions 


Three approximations have been made in the 
derivation of initial velocity steady state expression, 
equation (A2.15), and each of these will be dealt with in 
turn - generally and then with a specific case. 
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In the first approximation where d[HRP-II]/dt 
is set equal to zero and equation (A2.09) is derived, 


it is implicitly assumed in equation (A2.11) that 


4- ; 
= A lentl 
(ad [{HRP-II])/(dt Rare eae ]y) << [HRP] Equivalently 
4- : 

: a = 
ce Ns pape SNE Jo Ko obs a pseudo first order rate 
constant, then it is assumed: 

a [HRP-IT ] 

eae ary (AZ.17) 

at Ko obs HERP I, 


Reactions (A2.01), described by the differential 
equations (A2.04-2.06), under pseudo first order conditions 
where each of the substrates, hydrogen peroxide and ferro- 
cyanide, are in large excess of HRP and HRP compounds, are 
Settnre for Goncentrations and concentration derivatives at 
any time. A more complete description of the procedure 
is given by Evett (1968)*. The relative concentrations are 
pot red iit igure WAZ. 01) vs: time for agtypical case at 
a pH of 4.79. The following are the pseudo first order 
rate constants used in the calculation of the time course of 
the reaction progress curves of Figure (A2.01): 


(io 6 M0) ueeeck Crome Tm) 


(A2.18) 


k k [H, 


fobs 7° \leapome 27 0. 


‘Le Gi 10° sec 


oe ane) 106° Mimecece-) (3 .hesctlG aan) 


Koobs = jap 


MEE G pc (A2.19) 


i 





* 
I am indebted to Dr. M. Evett for his assistance in supply- 
ing a computer solution for the evaluation of the time 
course of the three component reaction system, 
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The apparent second order rate constants Kapp and arg 
were obtained from data collected in Chapter 3. The value 
of Kapp was obtained from Dolman (1969) which corresponds 
elVeceliy, tO Veluccmoptained byiChance (1963a).,.In Table 
(A2.01) are also listed normalized concentration deriva- 
tives at a number of representative reaction times for the 
reaction courses described in Figure (A2.01). 

Assume then that equation (A2.17) is valid when 


eoly some stialbl fraction, F,*of B.ts attained ‘or: 


d [HRP-TI | 


- = F (Ae) 
ele Rt ee 
obs a eo 


For 99% attainment of the steady state condition for this 


approximation or when F = 0.01: 
a [ERP= br 2 =| 
aeTHEP vs at (OOM ge “Sine eae) Sec 
= (P22) 
=45.12 SEC 


From Table (A2.01) a value of 3.2 sec? for a@[HRP-II]/(dt 
[HRP] 9) is yveached in eapproximately JS msec. “Similarky at 
can be shown that equation (A2.10) in equation (A2.11) 
involves the approximation: 


d [HRP ] 


= F (A2.23) 
SR Tae 
dt KI ops LHRP 0 


and if again Fe= 0.01 


fmeeeeO® Mo Mseo 7) (4.6) xet0r mM) 
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Table A2.01 


Normalized Concentration Derivatives Under 
Pseudo First Order Conditions at pH 4.79 


for Reactions (A2.01) 


: a [HRP] Ola Re aid) a{[HRP-IT] 
ae dt [HRP], dt THRPT ) sre aoe 
(scent) (see uw (sec *) 
ih acaba eer 107 Pea iee Vai 
2 Sais eeio. Apes SOG - Hag BOTH? 
4 Sey se ke ie oo. fail se Ne 
6 2 oe alee SS TO Soy ee 
8 Oe non -3.0 x 10+ 3.0 x lot 
10 au Soa ays Sonar Hat Seen 
12 9) Seen ma oe ae oy Se aa 
14 ase 107 - =f Pea 4.2 xtc? 
16 0 Se ye pl ee ae” 
18 aoe 108" ug hit cere Eevee kn 
20 si CO es a ity ii Bath so tee 
30 Gn Oa BENS abi Deal Oia 


40 DEO mn Oe =, eset GeO. sn Oi 


cnt 






ieooen.) bos tlaendyvouw 
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— — a -_ ee 


= . 
d LO7 = 
o _— - —_ SS + 






7" snlJ 
(seem) 


da3 


G[HRP]  _ 3 251 
dt THRPT , marty a hie | Oot LO sce 


25) (A2.24) 
lonsec 


From Table (A2.01), HRP reaches 99% of steady state con- 
ditions when the reaction has preceeded for only 3 msec, an 
even shorter time than for HRP-II as previously shown. 

The steady state approximation is used a third 
ender nat cime ena that as “in equation’ (A2. 13) "which is 


derived from equation (A2.12) by assuming that: 





3- 
d [Fe (CN) ] 
a[HRP-IT] 6 
aa aa << ay eee (AZE 25) 
Of 1£ O(HRP-Il]/dt has obtained some fractional portion, F, 
d {Fe (CN) 27) 
oem ars Ea the result may be written as: 
Era eeveny a 
dt [HRP], dt [HRP], j 
or: 
GQ[HRP-II] _ Fv 
dt [HRP], F [HRP], (A2.27) 


The value of V/THRP] 9 is simply the slope cfiapplomoiay 
vs. [HRP], aurea pi Of74 279. 1n Chapterpejs (Figures. 0S )meane 
has been determined to be 42 sec + by a weighted linear 
least squares analysis. If F = 0.01 again or the steady 
state conditions for this approximation are 99% attained 
then: 


qd [HRP-IT} 1 


= 0.01 x 42 sec 
ah ee 
GEIHRP) 5 (A2.28) 
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a{HRP-IT } 


BPEOMG Tab lem(henO eens valine: of dt (HRPT 5 


is attained at 
approximately 21 msec after the start of the reaction. The 
last approximation takes the longest to become valid and 

thus determines the attainment of the steady state conditions- 
99% attained 21 msec after the start of the reaction. The 
value of 21 msec is a much shorter time than that over which 
the initial velocities were calculated (a few seconds) and 


hence the steady state approximations are valid for the 


determination of v. 
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Appendix 3 
Purification of HRP by Sephadex Column Chromatography 


A procedure was developed for the purification 
of HRP in Sephadex column chromatography. Sephadex is a 
high molecular weight polysaccharide that is highly cross- 
linked. Molecules larger than’ the pores present in 
swollen Sephadex do not penetrate and pass through the 
column first. Smaller molecules penetrate into the pores 
and are eluted last. Molecules, in general, are eluted 
from a Sephadex column in order of decreasing molecular 
weight. 

The Sephadex used in this purification of HRP 
was a cation exchange gel, Sephadex C-50 (CM), and contains 
aelaLlge number Of (Carboxylucracidagroups. tsgAtna pHrvalue 
larger than the pK's of the carboxylic acid groups the acid 
groups are negatively charged and absorb positively charged 
solutes. Thus, separation of solute components on Sephadex 
C-50 (CM) is due partly to charge effects and to the mole- 


cular weight of the solute. 


Experimental Procedure 

Sephadex C-50 (CM) (7.5 gm) was equilibrated for 
48 hours in? 0.1 tonic strength phosphate buffer of pH 6.2 
with several changes of supernatant to ensure complete 
equilibration with buffer. The column (Sephadex Laboratory 


Column type K25/45; 2.5 x 45 cm) was equipped with a gradient 
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apparatus to give a continuously linearly increasing ionic 
strength gradient. With the gradient apparatus the ionic 
strength increases linearly from,0.1 to 0.2. 

The HRP (purified, lyophilized from Boehringer- 
Mannheim) (0.5 gm) was dissolved in 5 ml of water and added 
to the top of the column and eluted with a pressure head of 
about 5 cm of water. About 60 ml of effluent were allowed 
to pass through the column before the first samples were 
collected. Generally the major peak was observed to pass 
through the column at an effluent volume of 90 ml. In 
addition to the major peak two small poorly separated peaks 
remained at the top of the column. These components did 
not pass through the column even after a large volume of 
0.2 ionic strength phosphate buffer had been applied. 

The effluent was monitored by taking successive 
5 ml samples and measuring their absorbance after a 10:1 
dilution at 403 my and 280 mu. The absorbance at 403 mu is 
due mainly to the heme prosthetic group of HRP and the 
absorbance at 280 mu is due mainly to the apoprotein. The 
ratio of the ‘absorbance at 403 my to that at 280 mu is called 
the R.Z. and is often taken as an index of HRP purity, the 
higher the R.Z. the greater the purity of the enzyme pre- 
paration. 

Absorbances at 403 and 280 mu as well as R.Z. 
GLesDLOL Led Reece re ol lient volun in igure (A3.01)), 
The elution profile is not flattened at the center indica- 


ting that the sample size was not excessive for the type 
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Figure A3.01: Absorbances at 403 and 280 mp and R.Z. 


plotted vs. effluent volume. 
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-of column used. 

In general the effluent volume samples from 75 
to 95 ml were combined together to give a sample of HRP 
of R.Z. 2.8 or higher. The combined samples were precipitated 
with an 80% saturated ammonium sulfate solution by spinning 
on a centrifuge at 4°C. The resulting ammonium sulfate pre- 
cipitate containing the HRP was stored under refrigeration 
Perma bOout ia. 

Before use the precipitate was dissolved an water 
(5 ml), dialyzed exhaustively against cold distilled water 


and centrifuged and filtered with a Millipore filter. 
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Appendix 4 


Derivation of the Steady State Relaxation Time for the 


Binding of Imidazole to Hemin 
For the reactions, 


M+ Im —22R. yim (A4.01) 
k 
= app 


k 
Nien os fe SUSE MIm, (A4.02) 


where M is the metalloporphyrin, ferriprotoporphyrin IX, 


Im is imidazole and MIm and MIm., are the monoimidazole and 


bisimidazole hemin complexes, all without regard to state 
of protonation of the reactants. The differential rate 


expression is: 


d[MIm] _ 
atc: Ky app Mf] [Im] + k_,[MIm,] - jhe 


(A4.03) 


- Espey EE aa 


With the steady state approximation d[MIm]/dt = 0 in 


equation (A4.03) and solving for [MIm]: 


Ve cee 8 Gh ed [MIm., ] 
iim] 2 aPP : qa : (A4.04) 


-lapp ti 2app 


The differential rate expression in terms of [MIm, ] 1S 3 


ad [MIm., ] 


Sar appr mt! tim] - ace. te) (A4.05) 





From reactions (A4.01) and (A4.02) it is also easily shown 


thats: 
2 Ky app*2app re [Mim, ] 


—— (A4.06) 
&PP E app seep [M] [Im]? 


where the barred concentrations are equilibrium concentra- 


erons. 


Assuming a perturbation from the equilibrium 
concentrations of the amount A[{Mim,], the conservation rela- 


tions are then: 


[M] = [M] - A[MIm, ] (A407) 
[Im] = [Im] - 2A [MIm, ] (A4.08) 
[MIm, ] = [MIm, ] + A[MIm,] | (A4.09) 


Substituting equation (A4.07-4.09) in equation (A4.05) 


d([MIm, ] = A[MIm, ]) 


dt Bi ona aoet Aim) Cm) — 


2A [MIm, 1) * x ({iiim,] + AfMtm]) } 


K_2app*-lapp 


{k + ee 
fi -lapp K 2 app ( im] 2 [MIm.]) } (A4.10) 
First making the approximation: 


[im] = [im] = 2\ [MIm, ] (A471) 
and 


(A [MIm, 1)? << A[MIm, ]} (A4.12) 
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and then substituting equation (A4.06) in equation (A4.10), 


the result. is: 


a A(MIm, ] k k [Im] z 


k k 
2app lapp ~Zapp -lapp (A4.13) 
dt : 


SSlae + Kapp t ms 
The reciprocal relaxation time is given by the identity 
(a (A [MIm,]) /at) /A(MIm,] = (1) 77 (A4.14) 


Thus with equations (A4.14) and (A4.06) equation (A4.13) 


reduces to: 


——.,2 
k 1+ K Tm 
app apo (A415) _ 


Alr 


Koa -—- 
Pep eeSkE rm] 
-lapp 
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